(2x  IIBBIS 

BBwnsnaaw 

wniM$ 


THE  UNIVERSITY  OF  ALBERTA 

RELEASE  FORM 

NAME  OF  AUTHOR  Shu-Ming  Albert  LAI 

TITLE  OF  THESIS  Monte  Carlo  Study  of  Pattern  Loadings  on 

Continuous  Beams  and  Slabs 

DEGREE  FOR  WHICH  THESIS  WAS  PRESENTED  Master  of  Science 

YEAR  THIS  DEGREE  GRANTED  Spring  1981 

Permission  is  hereby  granted  to  THE  UNIVERSITY  OF 
ALBERTA  LIBRARY  to  reproduce  single  copies  of  this 
thesis  and  to  lend  or  sell  such  copies  for  private, 
scholarly  or  scientific  research  purposes  only. 

The  author  reserves  other  publication  rights,  and 
neither  the  thesis  nor  extensive  extracts  from  it  may 
be  printed  or  otherwise  reproduced  without  the  author's 
written  permission. 

( 

P 


DATED  .  T2;  r^r 


1980 


THE  UNIVERSITY  OF  ALBERTA 


Monte  Carlo  Study  of  Pattern  Loadings  on  Continuous  Beams 

and  Slabs 


© 


by 

Shu-Ming  Albert  LAI 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  Master  of  Science 

Department  of  Civil  Engineering 


EDMONTON,  ALBERTA 
Spring  1981 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


https://archive.org/details/Lai1981 


THE  UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 


The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  and  Research, 
for  acceptance,  a  thesis  entitled  Monte  Carlo  Study  of 
Pattern  Loadings  on  Continuous  Beams  and  Slabs  submitted  by 
Shu-Ming  Albert  LAI  in  partial  fulfilment  of  the 
requirements  for  the  degree  of  Master  of  Science. 


Abstract 


A  design  load  pattern  is  proposed  to  represent  the 
0.999  fractiles  of  the  maximum  lifetime  live  load  effects. 
The  proposed  load  pattern  is  different  from  the  traditional 
ACI  load  pattern,  which  produces  maximum  load  effects  with 
inconsistent  probability  of  being  exceeded.  Monte  Carlo 
simulations  have  been  used  to  generate  families  of  the 
lifetime  maximum  loading  effects  in  a  series  of  continuous 
beam-column  frames.  Linear  programming  is  then  used  to 
compute  the  magnitudes  of  the  pattern  loads  for  design.  The 
effect  of  pattern  loadings  for  the  dead  load  is  found  to  be 
insigni f icant . 
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1  .  INTRODUCTION 


1.1  Pattern  loadings 

1.1.1  Pattern  Loadings  of  Live  Load 

In  structural  design,  it  is  not  the  design  load  but  the 
effect  of  the  load  (moment,  axial  force,  etc),  which  is  of 
ultimate  interest.  Design  codes  in  the  United  States  (1)  and 
Canada  (19)  require  that  structural  elements  shall  be 
designed  for  the  maximum  effect  of  the  design  live  load.  The 
distribution  of  loads  causing  the  maximum  load  effect  in  a 
particular  structure  can  be  determined  by  drawing  an 
influence  line  for  the  structure  and  applying  loads  on  those 
parts  of  the  structure  where  their  effects  will  be  additive. 
In  a  continuous  beam  or  a  one  way  slab,  the  influence  lines 
of  all  major  load  effects  (8)  indicate  that  the  loadings 
required  for  maxima  can  be  reduced  to  two  types: 

a.  full  loading  on  alternate  spans, 

b.  full  loading  on  two  adjacent  spans  and  on  alternate 
spans  beyond  these. 

These  two  are  used  interchangeably  as  required  to  give  the 
maximum  load  effect  at  each  section. 

Section  8.9  of  the  AC  I  code  (1)  follows  a  similar  load 
pattern  for  the  live  load  effect  except  that  the  second  type 
is  simplified  to  full  loading  on  two  adjacent  spans  only. 

The  change  has  little  effect  on  the  cor respondi ng  maximum 
load  effect  because  the  influence  ordinates  in  the  alternate 
spans  are  relatively  very  small.  The  ACI  load  pattern  is 
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referred  to  in  this  report  as  the  "traditional  load 
pattern" . 

The  probability  of  occurrence  of  the  above  extreme 
loading  cases  is  not  Known,  however.  Moreover,  this 
traditional  pattern  may  not  produce  the  maximum  load  effects 
with  a  consistent  probability  of  being  exceeded. 

1.1.2  Pattern  Loadings  of  Dead  Load 

The  common  practice  of  structural  design  in  the  United 
States  and  Canada  does  not  take  into  account  the  pattern 
loading  effects  of  the  dead  load.  In  other  words,  the 
factored  dead  load  is  considered  constant  on  all  spans  in  a 
structure.  Section  2.3.3. 1  of  the  British  concrete  design 
code  (4),  however,  requires  that  a  pattern  loading  of  dead 
loads  should  be  considered  with  the  minimum  design  load 
equal  to  the  nominal  dead  load  and  the  maximum  equal  to  the 
factored  dead  load  to  account  for  the  most  unfavorable 
condition.  The  load  pattern  considered  is  similar  to  the 
traditional  load  pattern  (Sect.  3.2.2. 1  of  Ref.  4). 


1.2  Purpose  and  Scope  of  this  Study 

The  purpose  of  this  study  is  to  find  a  simple  design 
live  load  pattern  which  will  produce  all  load  effects  of 
interest  with  a  consistent  probability  of  being  exceeded  for 
the  extreme  cases.  As  a  part  of  this  study,  the  traditional 
load  pattern  is  evaluated  statistically.  Simple  analyses  of 
a  series  of  elastic  frames  are  carried  out  to  compute  the 
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load  effects  correponding  to  randomly  generated 
distributions  of  live  loads.  A  Monte  Carlo  process  is  used 
to  generate  the  loadings  according  to  the  probabilistic  live 
load  model,  which  has  been  described  in  the  literature  (5, 

6,  11,  16,  20,  21).  Linear  programming  is  used  to  develop  an 
equivalent  design  loading  pattern  for  the  extreme  loading 
cases.  Although  the  Monte  Carlo  study  is  based  on  office 
live  loads,  the  methodology  can  be  extrapolated  to  other 
types  of  occupancy. 

For  the  dead  load,  the  purpose  of  the  study  is  to 
investigate  the  effect  of  pattern  loadings.  Because  of  the 
simplicity  of  the  probabilistic  dead  load  model,  direct 
statistical  calculation  of  dead  load  effects  for  the  extreme 
cases  is  performed.  Again, linear  programming  is  used  to 
develop  an  equivalent  design  dead  load  pattern,  which  is 
compared  with  North  American  and  British  practice. 


1.3  Outline  of  Contents 

Chapter  2  contains  statistical  descriptions  of  the  live 
load  and  the  dead  load,  which  have  been  developed  by  several 
investigators  (11,16,20).  The  analytical  study  is  described 
in  Chapter  3,  which  describes  the  Monte  Carlo  analysis  and 
linear  programming  for  the  live  load  effects,  and  the 
statistical  calculations  of  the  dead  load  effects.  The 
details  of  the  structures  studied  are  also  presented.  In 
Chapter  4,  the  results  of  the  study  are  discussed  together 
with  the  development  of  a  most  desirable  design  load  pattern 
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for  the  live  load  and  the  dead  load.  Finally,  a  summary  and 
conclusions  are  presented  in  Chapter  5. 


2.  LOADING  MODELS 


2.1  Live  Load  Mode 1 

2.1.1  Nature  of  Live  Load 

The  gravity  live  load  on  a  floor  area  can  be 
represented  by  the  two  components:  a  sustained  load  and  a 
transient  load.  The  sustained  load  acts  continuously  in  time 
and  remains  relatively  constant  between  distinct  changes 
which  are  generally  considered  to  occur  when  there  are 
changes  in  tenants  or  occupancy.  The  sustained  load  consists 
of  furniture,  moveable  property  and  personnel  normally 
present.  Such  loads  are  usually  measured  in  live  load 
surveys.  For  certain  periods,  referred  to  in  this  report  as 
"vacant  periods",  the  sustained  load  may  be  entirely  absent. 

The  transient  load  happens  infrequently  with  a 
relatively  high  intensity  and  short  duration  in  the  order  of 
hours.  It  is  caused  by  unusual  events  such  as  crowding  of 
people  and  stacking  of  furniture  on  a  certain  area.  These 
special  loading  situations  are  not  usually  observed  in  live 
load  surveys. 

2.1.2  Live  Load  Surveys 

A  major  live  load  survey  (18)  comprising  thirty  two 
buildings  was  carried  out  in  England  in  1965  to  1967  by  the 
Building  Research  Station  (BRS).  More  recently  an  extensive 
live  load  survey  (9)  of  office  buildings  in  the  United 
States  of  America  was  conducted  in  1974  to  1975  by  the 
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National  Bureau  of  Standards  (NBS).  Twenty  three  buildings 
were  surveyed.  Statistical  analysis  (9)  of  the  NBS  survey 
data  indicates  that  the  mean  sustained  live  load  is  strongly 
correlated  to  room  use,  but  is  independent  of  room  area.  In 
order  to  enhance  the  flexibility  of  room  usage,  the  mean 
load  of  all  office  buildings  surveyed  is  used  in  the  live 
load  model  presented  in  Section  2.1.3.  The  mean  unit  load  of 
all  the  randomly  selected  buildings  from  the  NBS  survey  is 
0.555  KN/square  meter  (11),  which  is  very  close  to  0.565 
KN/square  meter  (17)  obtained  from  the  BRS  survey. 

The  variance  of  unit  load,  Var(u),  appears  to  decrease 
with  increasing  floor  area,  A  (see  Table  7  of  Ref.  18, 

Figs  29,30  of  Ref.  9  and  Fig.  2  of  Ref.  11).  The  analysis 
(11)  of  the  NBS  survey  data  yields  the  following  equation, 
discussed  more  fully  in  Section  2.1.3(a): 

Var (u)  =  0.0601  +  1  .407/A  (2.1) 

In  comparison,  Var(u)  obtained  by  McGuire  and  Cornell  (17) 
from  the  BRS  survey  data  is: 

Var (u)  =  0.0466  +  1.782/A  (2.2) 

The  values  of  Var(u)  computed  from  these  equations  have  been 
shown  to  be  very  close  (see  Fig.  3  of  Ref.  11). 

The  results  of  other  surveys  associated  with  a  room 
area  of  about  18.6  square  meters  (200  square  feet)  are 
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listed  in  Table  2.1(5).  The  survey  area  weighted  average  of 
the  mean  unit  load  is  0.566  KN/square  meter.  Values  ranging 
from  0.555  to  0.575  have  been  presented  in  Refs.  11,  12  and 
17.  The  survey  area  weighted  average  of  the  standard 
deviation  based  on  an  area  of  18.6  square  meters  is  0.35 
KN/ square  meter . 

2.1.3  Sustained  Live  Load  Model 

(a)  Sustained  Load 

The  instantaneous  sustained  live  load  intensity, 
w(x,y),  at  any  point  in  a  building  can  be  represented  by  a 
1 i near  relation  (21): 

w(x,y)=m+'tfb+tff  +  € ( x , y )  (2.3) 

In  this  relation,  m  is  the  mean  survey  load  for  the  type  of 
occupancy  considered.  Two  independent  zero-mean  random 
variables,  and  .represent  the  deviation  of  the 
building  average  unit  load  from  m  and  the  deviation  of  the 
floor  average  unit  load  from  the  building  average,  m  +  Xb  , 
respectively.  The  zero  mean  property  of  ^b  and  is 
justified  as  the  mean  unit  load  is  independent  of  the  floor 
area  (Sect.  2.1.2).  The  term  £(x,y)  is  a  stochastic  process 
representing  the  deviation  from  the  floor  average,  m  + 

+  .  It  has  a  zero  mean  and  is  independent  of  X  terms. 

In  general,  the  values  of  6(x,y)  at  two  points  are 
correlated  because  if  the  load  intensity  is  higher  than  the 
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TABLE  2.1  SUMMARY  OF  LIVE  LOAD  SURVEY  RESULTS 
FOR  OFFICES  (From  Ref.  5) 


Surveyed 

area 

Survey 

load 

Year 

Survey 

PI  ace 

Mean 

Std.  dev. 

(m2) 

(KN/m2) 

(KN/m2) 

1893 

Blackal 1 

U.S.A. 

7100 

0.780 

- 

1923 

Col  ey 

U.S.A. 

3700 

0.555 

- 

B1 ackal 1 

U.S.A. 

1100 

0.302 

0.13 

McIntyre 

U.S.A. 

3470 

0.436 

- 

1931 

Whi  te 

Engl  and 

1  3300 

0.517 

- 

1947 

Dunham 

U.S.A. 

46  370 

0.685 

- 

1952 

Dunham 

U.S.A. 

1270 

0.484 

- 

et  al 

1968 

Bryson  & 

U.S.A. 

4990 

0.570 

0.25 

Gross 

4500 

0.464 

0.21 

1969 

Karman 

Hungary 

1  3900 

0.603 

0.33 

1970 

Mitchell  & 

Engl  and 

160,000 

0.565 

** 

0.38 

Woodgate 

1973 

Paloheimo 

Fi  nl  and 

3000 

0.330 

- 

1974 

Cul  ver 

U.S.A. 

53900 

0.555 

*** 

0.37 

Dayeh 

Austral i a 

28000 

0.412 

0.25 

Schwartz 

U.S.A. 

400 

0.776 

0.52 

Area-wei ghted 
avera  ge 

0.566 

0.35 

2 

*  Based  on  an  office  area  of  approximately  18.6  m  (200  sq.ft) 
**  Derived  from  Eq.  2.1 

***  Derived  from  Eq.  2.2 
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floor  average  at  a  certain  point,  then  it  is  likely  that  the 
load  at  a  nearby  point  is  also  high.  Nevertheless,  Hauser 
(13)  has  shown  that  this  correlation  decays  rapidly  as  the 
distance  between  the  loads  increases.  The  analysis  was  based 
on  the  BRS  survey  results  (18);  the  correlation  coefficient 
was  shown  to  be  less  than  0.25  for  distances  greater  than  3 
meters.  Therefore,  if  the  floor  area  is  not  too  small, 

€(x,y)  can  be  reasonably  assumed  to  be  an  uncorrelated 
process  (17). 

Based  on  Eq .  2.3  and  the  preceding  assumptions,  the 
mean  and  standard  deviation  of  the  i nstantaneuos  unit  load, 
u,  on  a  floor  area,  A,  are  (17): 

E ( u )  =  m  (2.4) 

Var ( u )  =  a  +  b/A  (2.5) 

in  which  a  is  the  sum  of  the  variance  of  \  and  Yf  ,  and  b  is 
an  experimental  constant. 

Eq .  2.5  suggests  that  when  A  is  sufficiently  large, 
Var(u)  approaches  a.  Hence  a  is  determined  from  survey  data 
for  a  large  area.  The  value  for  b  is  selected  such  that 
Eq.  2.5  gives  a  good  fit  to  survey  data  for  other  areas.  The 

BRS  survey  of  office  live  loads  shows  a  standard  deviation 

of  0.216  kN/square  meter  at  A  =  192  square  meters  for  floors 
other  than  basements  and  grounds  (Table  7  of  Ref.  18).  This 
yields  a  =  0.0466  kN2/m4 .  In  Ref.  17  the  BRS  data,  which  are 
presented  in  Table  7  of  Ref.  18,  are  approximated  with  a 
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2  ? 

curve  based  on  Eq.  2.5  with  b  =  1.782  KN  /m  .  The  analysis 
(11)  of  the  NBS  survey  data  of  office  live  loads  presented 
in  Figs.  29,  30  of  Ref.  9  yields  a  =  0.0601  KN2/m4  and  b 
=  1.407  KN  /m  .  The  two  curves  are  very  similar  as  shown  in 
Fig.  3  of  Ref.  11.  Using  the  standard  deviation  of  0.35 
KN/square  meter  for  A  =  18.6  square  meters  presented  in 
Section  2.1.2  and  setting  a  equal  to  3.0  %  of  b,  compared  to 
2.6%  from  the  BRS  survey  and  4.3%  from  the  NBS  survey, 
yields  a  =  0.044  KN2/m4  and  b  =  1.460  kN2/m2.  These  values 
are  used  for  analysis. 

(b)  Sustained  Load  Effect 

To  relate  a  random  set  of  loads  distributed  over  an 
area  to  a  particular  load  effect  (moment,  shear,  etc.),  the 
sum  of  the  products  of  the  magnitude  of  the  loads  at  each 
point  times  the  height  of  the  influence  surface  for  that 
load  effect  at  that  point  is  used.  The  mean  and  standard 
deviation  of  the  sustained  equivalent  uniformly  distributed 
load  (EUDL),  L,  which  will  produce  the  same  load  effect  as 
the  actual  random  set  of  loads  are  (17): 


E  (  L  )  = 

m 

(2.6) 

Var ( L ) 

=  a  +  Kb/ A  j 

(2.7a) 

in  which 

k= 

'o 

(VuWff 

'0  '  <  o 

r 1  I2 

I(x,y) 
o  ; 

(2.7b) 

The  function 

I ( x , y )  is  the 

norma  1 i zed 

influence  surface  of 

the  load  effect  of  interest  over  an  influence  area,  A(  . 

The  influence  area,  A(  ,  is  generally  assumed  to  be 
equal  to  the  floor  area  over  which  the  influence  surface  of 
the  load  effect  considered  is  significantly  different  from 
zero  (17).  For  beam  load  effects  in  a  conventional  one  way 
slab  or  beam-column  frame  building,  the  significant 
influence  area  is  generally  considered  as  twice  the 
conventional  tributary  area  of  one  span.  For  the  column  load 
effects,  the  significant  influence  area  is  4  times  the 
tributary  area. 

The  values  of  k  have  been  shown  to  be  insensitive  to 
different  end  conditions  of  a  single  beam  and  number  of 
spans  in  a  frame.  Using  approximate  polynomial  influence 
surfaces,  k  is  found  to  be  2.04  for  end  moments,  2.76  for 
mid-span  moments,  2.20  for  column  loads  and  similar  values 
for  other  load  effects  (17).  Thus  the  variance  Var(L) 
computed  using  Eq .  2.7(a)  for  A!  =37.2  square  meters  (400 
sq .  ft.)  ranges  from  0.144  to  0.179  for  the  extreme  cases. 
Because  of  the  relative  insensitivity  of  k  to  load  effect 
types,  k  =  2.2  for  all  load  effect  types  has  been  suggested 
(11)  and  it  is  used  in  this  study. 

(c)  Statistical  Distribution  of  Sustained  Load 

Five  major  live  load  surveys  have  been  analysed  in 
Ref.  6  using  three  probability  models:  normal,  lognormal, 
and  gamma.  The  results  indicate  that  the  gamma  distribution 
is  generally  better  in  the  overall  fit  and  gives  the  best 
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fit  in  two-thirds  of  the  cases.  Furthermore,  the  gamma  model 
is  in  good  agreement  with  the  survey  data  of  most  occupancy 
types  in  upper  tail  region  and  even  at  the  99.9%  level  for 
offices  (see  Table  2  of  Ref.  6).  As  the  EUDL  differs  from 
the  actual  unit  load  only  by  the  weighting  function  I(x,y), 
the  instantaneous  EUDL  is  also  considered  as  gamma 
di str i buted . 

(d)  Statistical  Distribution  of  Sustained  Load  Duration 

Throughout  the  life  of  a  floor  area,  the  sustained  load 
is  likely  to  be  influenced  primarily  by  the  change  in 
occupancy  while  minor  fluctuations  during  the  same  occupancy 
can  be  ignored.  Therefore  the  sustained  load  is  assumed  to 
be  constant  during  the  time  between  occupancy  changes.  The 
load  during  any  such  constant  period  will  be  represented  by 
the  instantaneous  sustained  EUDL. 

The  NBS  office  live  load  survey  data  for  occupancy 
duration  exhibit  an  exponential  distribution  as  shown  in 
Fig.  18  of  Ref.  9.  The  mean  duration  is  8.0  years,  which  is 
close  to  the  mean  value  of  8.8  years  obtained  by  the  BRS 
survey  (18).  The  BRS  survey  data  for  occupancy  duration  also 
agree  fairly  well  with  an  exponential  model  as  shown  in 
Fig.  5  of  Ref.  21.  Figure  18  of  Ref.  9  suggests  a  minimum 
duration  of  1  year.  In  this  study  the  the  duration  of 
occupancies  is  modelled  with  an  exponential  distribution 
with  a  minimum  duration  of  1  year  and  a  mean  of  8  years. 
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2.1.4  Vacant  Period 

Vacancy  or  unloading  events  in  an  area  are  generally 
neglected  in  deriving  the  lifetime  maximum  EUDL  as  a 
conservative  assumption  (17).  In  this  study,  however, 
unloading  events  in  a  span  are  significant  and  should  be 
considered.  Although  no  data  are  available  about  office 
vacancy,  it  is  likely  that  unloading  occurs  during  periods 
of  occupancy  changes,  repainting  or  remodelling  of 
furniture.  In  this  study,  an  area  is  arbitrarily  assumed  to 
be  entirely  unloaded  during  each  change  of  occupancy.  Within 
the  same  occupancy,  additional  unloading  events  are  assumed 
to  occur  every  2  years. 

The  duration  of  an  unloading  event  is  assumed  from  1 
day  to  2  months  and  is  approximated  by  a  uniform 
distribution.  The  sensitivity  of  load  effects  to  the 
above-mentioned  parameters  is  examined  in  Sect.  4.1.1(b). 

2.1.5  Transient  Live  Load  Model 

(a)  Transient  Load 

Very  few  data  are  available  on  the  transient  live  load. 
Nevertheless,  a  model  has  been  proposed  (20)  assuming  a 
series  of  randomly  distributed  load  cells,  each  of  which 
contains  a  cluster  of  concentrated  loads.  Similar  to  the 
sustained  load,  the  load  effect  is  also  considered  to  obtain 
the  EUDL.  The  mean,  mE ,  and  variance,^  ,  of  the  EUDL,  E 
associated  with  one  transient  load  event  are  (17): 


) 
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mE  =  AmpjmQ/Aj 

=AK(mRq-2  +  h1qTr  + 


(2.8) 

(2.9) 


in  which  A=  mean  number  of  load  cells  in  the  influence 
area,  A(  .  Q  and  R  are  random  variables  representing  the 
weight  of  of  a  single  concentrated  load  in  the  cell  and  the 
number  of  loads  per  cell  respectvely.  K  is  equal  to  2.2 
(Sect.  2.1.3(b)). 

There  are  no  data  available  with  which  Eqs .  2.8,  2.9 
can  be  estimated  with  confidence.  McGuire  and  Cornell  (17) 
have  estimated  the  values  of  the  following  parameters 
(Table  1  of  Ref.  17):  (mQ,  ?Q  )  =  (0.65,  0 .  1 3 )  KN ;  (mR,^TR  ) 

=  (5,  2).  In  estimating  A,  which  is  area-dependent,  it  is 
reasonable  to  assume  that  the  average  number  of  load  cells 
per  unit  area  (and  therefore  the  mean  and  standard  deviation 
of  unit  load ).  decreases  with  increasing  area.  McGuire  and 
Cornell  (17)  have  proposed  the  following  values  of  A,  which 
are  a  function  of  influence  area,  A(  : 


Al 

>  37.2  m2 

A  =  7(4,  -  1 5 . 3  )  /  0 . 84 

A, 

=  27.9  m2 

A  =  3.99 

(2.10) 

Al 

=  18.6  m2 

A  =  2.7 

A, 

=  9.3  m2 

A  =  1.4 

The  mean  and  standard  deviation  of  the  unit  transient  load 
from  Eqs.  2.8  and  2.9  based  on  McGuire  and  Cornell's 
assumptions  are  plotted  as  a  function  of  influence  area  in 
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Fig.  2.1.  The  figure  shows  the  decrease  of  the  mean  and 
standard  deviation  with  increasing  area. 

Ellingwood  and  Culver  (11)  estimated  that  (m  ,  ^TQ  )  = 
(0.67,  0. 1 1 )KN,  (mR,  <JR  )  =  (4,  2)  and 

A  =  y ( A,  -  14.41/0.58  (2.11) 

Values  of  the  mean  and  standard  deviation  from  Eqs .  2.8  and 
2.9  based  on  these  assumptions  are  plotted  in  Fig.  2.1  for 
comparison  with  McGuire  and  Cornell's  curves.  The  figure 

shows  that  the  two  sets  of  assumptions  give  curves  that  are 

2  2 

very  similar  for  A!  >  30  m  .  For  A!  <  30  m  ,  Ellingwood's 
curve  appears  unreasonable  as  the  mean  unit  load  decreases 
very  rapidly  with  decreasing  area  (the  same  unreasonable 
trend  occurs  in  standard  deviation  as  A!  <  20  m  )  In 
personal  communication  with  Ellingwood,  it  has  been 
suggested  that  the  mean  unit  load,  and  therefore  the  average 
number  of  load  cells  per  unit  area,  should  be  assumed  to  be 
constant  when  the  influence  area  is  less  than  30  m2 ,  ie, 

A,  <  30  m2  A  =  0.  173A,  (2.12) 

The  mean  and  standard  deviation  based  on  Eq  2.12  are  also 
plotted  in  Fig.  2.1,  which  shows  a  reasonable  trend  in  the 
standard  deviation.  In  this  study,  Eqs.  2.11  (for  A!  >  30  m2 ) 
and  2.12  are  used  to  calculate  the  mean  and  standard 
deviation  of  the  transient  EUDL. 


Standard  Deviation  of  Transient  Load  Mean  of  Transient  Load 

KN/m2  KN/m2 
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Fig. 2.1  Mean  and  Standard  Deviation  of  Transient  Load 
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(b)  Statistical  Distribution  of  Transient  Load  Event 

Peir  (20)  has  shown  by  means  of  a  numerical  example 
that  the  distribution  of  the  transient  load  can  be 
represented  by  a  gamma  distribution.  The  occurrence  of  the 
transient  load  is  generally  assumed  to  follow  a  Poisson 
process  (17)  based  on  the  assumption  that  transient  load 
events  are  independent  in  space  and  time.  The  time  between 
two  consecutive  transient  load  events  is  assumed  to  conform 
to  an  exponential  distribution  with  a  mean  of  1  year  (11,17) 
for  offices. 

The  duration  of  each  transient  load  event  is  assumed  to 
be  8  hours  (15).  Moreover,  transient  load  events  are  assumed 
not  to  overlap,  although  a  transient  load  event  can  be 
followed  immediately  by  another  transient  load  event  so  that 
the  minimum  time  between  occurrences  is  8  hours. 

2.1.6  Live  Load  Model  Used  in  this  Study  -  Summary 

The  sustained  EUDL  is  gamma-distributed  with  a  mean  of 
0.566  KN/square  meter  and  a  variance  equal  to 

Var (L)  =  0.044  +  3.212/A,  (2.13) 

The  duration  of  a  sustained  load  event  is  modelled  with  an 
exponential  distribution  with  a  minimum  duration  of  1  year 
and  a  mean  of  8  years. 

One  unloading  event  occurs  at  the  start  of  each 
sustained  load  event  and  additional  unloading  events  occur 
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every  2  years  during  that  sustained  load  event.  The  duration 
of  a  vacant  period  is  approximated  by  a  uniform  distribution 
with  a  minimum  duration  of  1  day  to  a  maximum  of  2  months. 

The  transient  EUDL  is  assumed  to  be  gamma-distributed 
with  a  mean  equal  to 


mE  =  2.68  7\/ Aj 

and  a  variance  equal  to 

STE2  =  1 9 . 9A/A2 
where 

A  =  *J  ( A,  -  14.41/0.58 
A  =  0.173  A, 


(2.14) 


(2.15) 

for 

Al 

> 

30  m2 

(2.16) 

for 

A, 

« 

30  m2 

(2.16) 

The  duration  of  each  transient  load  event  is  8  hours.  The 
time  between  two  consecutive  transient  load  events  is 
modelled  with  an  exponential  distribution  with  a  minimum 
time  of  8  hours  and  a  mean  of  1  year. 


2.2  Dead  Load  Model 

The  dead  load  includes  the  weight  of  structural 
members,  permanent  equipment  and  permanently  supported 
non-structura 1  elements  such  as  partitions,  roofing  and 
installations.  During  the  life  of  a  structure,  the  dead  load 
maintains  a  relatively  constant  magnitude  though  slight 
changes  may  occur . 
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Magnitudes  of  dead  loads  may  vary  from  those  assumed 
due  to  variations  in  size  of  members,  density  of  material 
and  weight  of  non-structura 1  items.  Generally,  the  weights 
of  non-structura 1  elements  have  a  strong  effect  on  the 
variability  in  the  dead  load  (12). 

The  probability  distribution  of  the  dead  load  is 
generally  assumed  to  be  normal.  The  ratio  of  mean  load  to 
nominal  dead  load  is  commonly  assumed  to  be  unity  and  the 
coefficient  of  variation  (12)  from  0.06  to  0.15.  In  this 
study,  the  mean  load  is  assumed  to  be  1.05  times  the  nominal 
load,  and  the  coefficient  of  variation  of  the  dead  load  is 
taken  as  0.10,  based  on  Ref.  12. 


3.  ANALYTICAL  STUDY 


3.1  Live  Load  Study 

A  Monte  Carlo  analysis  is  used  to  generate 
distributions  of  the  maximum  lifetime  loadings  at  various 
points  in  a  structure.  A  linear  programming  technique  is 
then  used  to  select  a  pattern  loading  to  represent  a 
particular  fractile  of  the  maximum  lifetime  loadings. 

3.1.1  Monte  Carlo  Analysis 

(a)  Description  of  Monte  Carlo  Technique 

The  Monte  Carlo  technique  is  a  method  to  generate  a 
large  number  of  hypothetical  samples  by  computer 
simulations.  This  method  requires  a  mathematical 
relationship  between  the  variable  being  studied  and  each 
basic  variable  which  affects  it,  plus  a  statistical 
description  of  each  basic  variable.  A  set  of  values  of  the 
basic  variables  is  randomly  selected  according  to  their 
statistical  properties  and  the  final  variable  is  determined 
using  the  relationship  between  the  basic  variables  and  the 
final  variable.  This  procedure  is  repeated  many  times  to  get 
a  large  sample  of  the  final  variable  for  statistical 
analysi s . 

In  this  study,  the  Monte  Carlo  technique  is  used  to 
generate  large  samples  of  lifetime  minimum  and  maximum  live 
load  effects  at  various  points  in  a  given  frame.  During  the 
lifetime  of  the  frame,  random  values  of  the  sustained  load, 
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vacant  period  and  transient  load  are  continuously  generated, 
and  load  effects  are  calculated  accordingly.  The  minimum  and 
maximum  values  of  each  load  effect  during  the  lifetime  are 
eventually  selected.  This  procedure  is  repeated  enough  times 
(500  to  5000  times)  to  give  usable  statistical  distributions 
of  the  lifetime  minimum  or  maximum  values. 

(b)  Assumptions  in  the  Computer  Program 

The  computer  program  studies  a  given  one-story  frame 
that  is  randomly  loaded  according  to  the  live  load  models 
described  in  Chapter  2.  The  life  of  the  frame  is  assumed  to 
be  50  years. 

During  the  lifetime  of  a  given  floor  area,  successive 
loading  events  are  assumed  to  be  independent.  Based  on  the 
insignificant  correlation  between  the  sustained  live  loads 
of  two  rooms  found  in  the  analysis  (7)  of  loading  data,  it 
is  reasonable  to  assume  no  correlation  for  the  sustained 
load  before  or  after  an  occupancy  change.  Consecutive 
transient  load  events  (Sect.  2.1.5(b))  are  also  assumed  to 
be  independent,  and  so  are  unloading  events. 

The  number  of  rooms  and  room  arrangement  on  the  floor 
area  of  a  frame  give  rise  to  different  loading  changes 
occur ing  over  the  floor.  It  has  been  shown  (16)  that  it  is 
conservative  to  assume  the  entire  'significant7  influence 
area  to  be  occupied  by  one  room.  As  mentioned  in  Section 
2.1.3(b),  the  significant  influence  area  is  the  area  over 
which  the  influence  surface  of  the  load  effect  considered  is 
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'significantly'  different  from  zero.  For  most  beam  load 
effects,  the  influence  surface  in  the  spans  other  than  the 
span  where  the  load  effect  is  considered  is  relatively 
small.  For  this  reason  McGuire  and  Cornell  (17)  have  taken 
the  significant  influence  area  for  all  beam  load  effects  to 
be  the  influence  area  of  the  span  considered  (ie  twice  the 
tributary  area  of  the  span  considered). 

For  the  interior  support  moments,  however,  the 
influence  surface  is  usually  significant  in  the  two  spans 
adjacent  to  the  support.  Conservatively,  the  significant 
influence  area  equal  to  twice  the  tributary  area  of  two 
adjacent  spans  should  also  be  investigated.  This  matches  the 
situation  in  which  one  room  occupies  the  entire  influence 
area  of  two  adjacent  spans  so  that  the  same  loads  occur  over 
two  spans  at  any  time.  For  axial  forces  in  interior  columns, 
the  significant  influence  area  is  also  equal  to  twice  the 
tributary  area  of  two  spans  (ie  4  times  the  tributary  area 
of  the  column)  because  the  influence  surface  is  significant 
in  two  adjacent  spans.  McGuire  and  Cornell  (17)  have 
suggested  this  influence  area  in  the  case  of  interior 
columns  but  not  in  the  case  of  negative  moments  at  interior 
supports . 

In  the  Monte  Carlo  analysis,  two  loading  distributions 
are  assumed: 

Case  MCI  -  the  influence  area  of  each  span  is  occupied 


by  one  room. 


23 


Case  MC2  -  similar  to  Case  MCI  except  that  the  influence 
area  of  two  particular  adjacent  spans  are  occupied 
by  one  room. 

Because  Case  MCI  is  a  conservative  assumption  for  most  load 
effects  as  discussed  above  and  appears  most  likely  to 
happen,  it  is  investigated  for  all  load  effects  of  interest. 
For  comparison  Case  MC2  is  also  examined  for  interior 
negative  support  moments  and  interior  column  loads. 

Case  MCI  implies  that  all  sustained  load,  transient 
load  and  unloading  events  on  adjacent  spans  are  independent 
(Sect.  2.1.5(b)),  and  the  influence  area  used  in  Eqs .  2.7, 
2.9  to  calculate  the  standard  deviation  of  load  magnitudes 
is  twice  the  tributary  area  of  the  loaded  span.  For 
Case  MC2 ,  the  loads  are  identical  on  two  adjacent  spans 
where  the  interior  support  moments  and  the  column  load  are 
considered;  the  influence  area  used  in  Eqs.  2.7,  2.9  for  the 
loads  on  the  two  spans  is  twice  the  tributary  area  of  the 
two  spans.  The  loadings  on  other  spans  are  independent  as  in 
Case  MCI . 

The  computer  program  based  on  Case  MCI  is  described  in 
the  next  section.  For  Case  MC2 ,  the  computer  program  is  only 
slightly  modified  to  apply  identical  loads  on  two  particular 
adjacent  spans.  The  pattern  loading  models  developed  later 
will  be  based  on  Case  MCI  loadings  for  all  load  effects.  For 
comparison,  the  models  will  also  be  based  on  Case  MC2 
loadings  for  load  effects  which  are  more  critical  in 
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Case  MC2  loadings  (interior  column  loads,  etc.)  and  based  on 
Case  MCI  loadings  for  all  other  load  effects.  The  latter 
case  is  later  referred  as  Case  MC3. 

(c)  Description  of  the  Computer  Program 

The  computer  program  is  capable  of  generating  estimates 
of  the  lifetime  minimum  and  maximum  values  of  a  number  of 
load  effects  in  a  population  of  500  to  5000  identical  single 
storey  frames.  The  50  year  life  of  the  frame  is  divided  into 
8-hour  intervals  selected  so  that  the  length  of  one  interval 
is  the  same  as  the  duration  of  a  transient  load  with  the 
result  that  during  any  interval  the  loads  on  any  span  do  not 
change.  The  loads  on  each  span  are  randomly  chosen  from  the 
loading  distributions  and  the  load  effects  caused  by  the 
loads  occuring  in  each  interval  are  determined.  When  this 
has  been  done  for  each  load  interval  in  the  50  year  life, 
minimum  and  maximum  values  for  each  load  effect  during  this 
period  are  selected  and  denoted  as  the  lifetime  minimum  and 
maximum  values,  respectively.  The  above  steps  are  denoted  as 
one  run  in  the  computer  program.  In  most  cases  studied,  5000 
runs  for  one  particular  frame  are  used  to  obtain  5000 
lifetime  minimum  and  maximum  values  of  each  load  effect 
considered . 

The  structural  analysis  in  the  program  is  simplified  by 
using  an  influence  coefficient  matrix.  The  matrix  for  a 
particular  frame  is  multiplied  by  the  vector  consisting  of 
values  of  the  uniformly  distributed  load  on  each  span  at  a 
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particular  instant  to  obtain  the  corresponding  load  effects. 
The  formulation  of  the  influence  coefficient  matrix  is  shown 
in  Appendix  A. 

Condensed  flow  diagrams  of  the  computing  program  are 
shown  in  Figs.  3.1  to  3.3  (the  program  listing  is  in 
Appendix  B).  All  the  required  data  such  as  the  statistical 
parameters,  the  influence  coefficient  matrix  for  the  frame 
studied,  the  number  of  spans  and  the  number  of  runs  are 
first  read  in.  Then  the  random  sustained  loads,  unloading 
events  and  transient  loads  are  generated  for  each  span  in 
the  50  year  lifetime.  For  efficiency,  actual  calculation  of 
the  total  load  on  each  span  is  done  only  when  the  load  on 
that  span  changes;  the  interval  number  at  the  time  of  change 
is  also  recorded.  Similarly  the  load  effects  are  calculated 
using  the  influence  coefficient  matrix  whenever  the  total 
load  changes  on  any  span.  At  the  final  stage  of  each  run, 
lifetime  minimum  and  maximum  values  of  load  effects  are 
selected.  Finally,  the  lifetime  minimum  and  maximum  values 
of  each  load  effect  from  each  run  are  sorted  and  printed 
out . 

3.1.2  Linear  Programming 

(a)  Description  of  Linear  Programming 

Linear  programming  is  a  technique  to  provide  optimum 
solutions  to  problems  with  a  number  of  possible  solutions. 
The  problem  has  to  be  formulated  into  a  linear  program  which 
is  a  mathematical  model  consisting  of  an  objective  function 
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FIG.  3.1  CONDENSED  FLOW  DIAGRAM  OF  THE  MONTE  CARLO  PROGRAM 


- 
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SUBROUTINE  SUSTAN 


FIG.  3.2  SUBROUTINE  SUSTAN  OF  THE  MONTE  CARLO  PROGRAM 
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SUBROUTINE  TRANST 


FIG.  3.3 


SUBROUTINE  TRANST  OF  THE  MONTE  CARLO  PROGRAM 
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and  constraints,  both  in  terms  of  linear  expressions  of 
variables.  The  objective  function  is  the  factor  to  be 
optimized  such  as  maximization  of  profits  or  minimization  of 
costs.  The  optimization  is  restricted  by  the  constraints  on 
the  values  of  the  variables.  The  constraints  are  in  linear 
forms  of  inequality  expressions  or  equations. 

A  linear  program  with  two  variables  can  be  solved 
graphically.  For  more  than  two  variables  the  Simplex  method 
(10)  has  been  developed.  The  method  is  an  iterative 
technique  performed  using  a  computer.  Theoretically 
speaking,  the  method  is  capable  of  giving  the  true  optimal 
solution  to  the  problem  solved.  In  practice,  however,  the 
solution  may  be  inaccurate  due  to  the  rounding  errors 
accumulated  in  repeated  computations.  For  these  and  other 
reasons,  the  so  called  Revised  Simplex  method  has  been 
developed  to  reduce  the  errors  and  thus  improve  the 
solution.  This  method  is  described  more  fully  in  Ref.  10.  In 
this  study  the  Revised  Simplex  method  is  used  to  solve 
linear  programs;  the  computer  program  for  this  method  is 
available  from  the  IMSL  library  (14). 

In  this  study  the  Monte  Carlo  results  for  the  extreme 
loading  cases  are  approximated  with  design  load  patterns. 
Linear  programming  is  used  to  find  the  load  factors 
corresponding  to  a  particular  load  pattern  to  best 
approximate  the  Monte  Carlo  results.  The  objective  is  to 
select  load  factors  to  minimize  the  differences  between  the 
factored  load  effects  from  a  particular  load  pattern  and  the 
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corresponding  load  effects  from  the  Monte  Carlo  results.  A 
series  of  different  load  pattern  schemes  are  considered  to 
select  the  most  desirable  model. 

(b)  Design  Loadings  and  Load  Fractiles  Used  in  the  Study 

In  any  design  load  pattern  to  be  discussed,  the  load  on 
each  span  is  the  product  of  a  load  factor  and  the  nominal 
live  load.  The  nominal  live  load,  Ln ,  is  a  defined  quantity 
based  on  Section  4.1.6  of  National  Building  Code  of 
Canada  (19): 


for  AT  >  20  m2 
for  AT  4  20  m2 


Ln  =  (0.3  +n/9.8/At)L0 
Ln  =  L0 


(3.1) 


in  which  AT  is  the  tributary  area  in  square  meters.  Since 
consideration  is  given  to  office  buildings  in  this  study,  L0 
is  equal  to  2.4  KN/square  meter. 

Basically  two  types  of  design  load  patterns  should  be 
considered.  The  first  type  is  similar  to  the  traditional 
load  pattern  discussed  in  Section  1.1.1,  except  that  there 
may  be  a  small  load  on  the  unloaded  spans.  The  second  type 
is  a  simpler  pattern  suggested  by  Beeby  (2),  consisting  of 
heavy  loading  on  only  one  span  and  light  loading  on  all 
other  spans.  Other  patterns  to  be  considered  will  be  derived 
from  these  two  basic  types.  The  design  load  patterns  will  be 
described  more  fully  in  Chapter  4. 

In  order  that  the  factored  load  pattern  selected  in 
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this  part  of  the  study  should  have  a  probability  of 
occurrence  equivalent  to  that  implied  by  the  National 
Building  Code  (19)  or  ANSI  (12)  load  factors  for  live  load 
on  a  simple  span,  the  probability  of  occurrence  of  various 
factored  live  loads  has  been  studied.  The  probabilty  that 
the  lifetime  maximum  load  will  be  less  than  or  equal  to  the 
factored  live  loads  based  on  current  NBC,  proposed  ANSI  and 
current  ANSI  are  listed  in  Table  3.1  for  three  different 
influence  areas.  The  values  in  the  table  have  been 
determined  by  using  a  Monte  Carlo  program  based  on  the  one 
described  in  Section  3.1.1,  except  that  the  number  of  spans 
has  been  set  equal  to  one  and  the  output  is  expressed  in 
terms  of  lifetime  maximum  values  of  the  live  load  instead  of 
load  effects.  The  results  indicate  that  the  factored  live 
load  values  are  close  to  the  0.999  fractile  except  for  small 
areas.  For  this  reason,  the  0.999  fractiles  of  the  lifetime 
minimum  or  maximum  load  effects  are  used  in  the  modelling  of 
pattern  loadings. 

(c)  Formulation  of  the  Linear  Program 

For  convenience  in  the  following  discussion,  Mpj  is 
used  to  represent  the  absolute  value  of  the  algebraic 
minimum  or  maximum  factored  load  effect  i  which  can  be 
obtained  at  a  critical  section  from  a  design  load  pattern, 
and  Mmj  is  the  cor respondi ng  absolute  value  of  a  high 
fractile  of  the  distribution  of  the  minimum  or  maximum  load 
effect  i  from  the  Monte  Carlo  analysis.  Note  that  the 
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TABLE  3.1  PROBABILITY  OF  THE  LIFETIME  MAXIMUM  LOAD  BEING 
LESS  THAN  OR  EQUAL  TO  THE  FACTORED  LIVE  LOAD 


AI 

* 

1-5Ln 

(NBC) 

★  ★ 

’■6Ln 

( Proposed  ANSI ) 

*** 

]-7Ln 

( Current  ANSI ) 

20  m2 

0.925 

0.956 

0.973 

50  m2 

0.9982 

0.9986 

0.9954 

100  m2 

0.9984 

0.9988 

0.9912 
to  0.9970 

*  Assume  Ay  =  0 . 5 Ay  in  Eq.  3.1 

**  Lp  =  (0.25  +^)Lo  for  Aj  >  400  sq.ft 

Ln  =  Lq  for  Aj  <.400  sq.ft 

D 

***  L  =  [  1  -  min  {0.0008Ay,  0.6,  0.23(1  +  — -)}]L 

o 

Assume  Ay  =  0 . 5 Ay 
D 

—  =  0.67  to  2.0 
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subscript  p  refers  to  pattern  and  m  refers  to  Monte  Carlo. 

In  the  optimization,  the  penalty  for  values  of  Mmi  less 
than  Mpi  (an  unconservative  case)  is  twice  that  for  Mmi 
greater  than  Mpj  .  The  objective  is  to  minimize  the  sum  of 
penalties  for  all  load  effects  to  be  fitted  by  a  load 
pattern.  In  this  way  the  results  tend  to  be  conservative. 

As  discussed  in  the  previous  section,  Mmj  are  taken 
equal  to  absolute  values  of  the  0.999  fractiles  from  the 
Monte  Carlo  results.  Two  cases  are  assumed  in  the  modelling 
with  a  load  pattern.  The  first  case  is  that  the  Monte  Carlo 
results  based  on  Case  MCI  discussed  in  Section  3.1.1(b)  are 
used  for  Mmi  in  the  linear  programming.  The  second  case  is 
that  the  Monte  Carlo  results  based  on  Case  MC3  discussed  in 
Section  3.1.1(b)  are  also  used  for  Mmj  . 

The  load  effects,  Mpj  ,  for  a  given  load  pattern  and 
frame  are  determined  in  the  linear  function  of  load  factors 
using  the  influence  coefficient  matrix  (Appendix  A)  as 
fol lows : 


Mpi  =  aLn  C-|  +  bLn  C2  (3.2) 

where  a  and  b  are  load  factors,  which  will  be  defined  more 
fully  in  Chapter  4.  The  nominal  live  load,  Ln ,  has  been 
defined  by  Eq.  3.1.  The  term  C1  is  the  sum  of  the  influence 
coefficients  of  load  effect  i  for  the  spans  which  are  loaded 
by  the  large  uniform  load,  aLn ,  and  C2  is  the  sum  of  the 
remaining  influence  coefficients  of  load  effect  i.  Because 
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MPj  are  defined  as  absolute  values  in  the  linear  program, 
the  algebraic  sign  of  each  influence  coefficient  for 
negative  load  effects  (negative  moments,  negative  shears) 
has  to  be  converted  to  the  opposite  sign  when  formulating 
Eq.  3.2. 

The  linear  program  is  formulated  as  follows:  The 
objective  function  is 

n 

minimize  z  =  ^>dj  (3.3) 

i=i 

where  n  is  the  total  number  of  load  effects  to  be  fitted, 


and  dj 

i  s  the 

penalty  for  load  effect  i. 

The  penal ty ,  d , , i s 

defined 

i  n  the 

fol lowing: 

If  Mpi 

>  Mm, 

di 

■  "pi 

i 

(3.4) 

If  Mpi 

Mmi 

di 

=  2(Mm 

i  -  Mpj  ) 

(3.5) 

Note  that  both  Mpi  and  Mmj  are  absolute  values.  The 
parameters  dj  and  Mpi  (in  terms  of  a  and  b)  are  variables, 
while  Mmi  are  Known  values  as  defined  previously.  Eqs.  3.4, 
3.5  are  transformed  into  constraints  in  the  linear  program 
as  shown  in  Eq .  3.6.  The  inequality  expressions  are  used 
instead  of  equations  because  eqs.  3.4  and  3.5  are  never 
satisfied  simultaneously.  This  penalty  expression  is 
arbitrarily  chosen  such  that  unconservative  estimates  of 
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load  effects  will  receive  twice  the  penalty  as  an  equal 
conservative  estimate.  No  significant  difference  in  load 
factors  has  resulted  from  using  multipliers  of  3  or  4  in 
Eq .  3.5. 

The  constraints  are: 

Mpi  ~  ^mi  ^  dj 

for  i  =  1  ,  n  (3.6) 

Mmi  ~  Mpj  ^  0 . 5dj 

The  variables,  a,  b  and  dj  ,  are  also  defined  as  non-negative 
variables  in  the  linear  program.  In  addition  to  the  above 
constraints,  the  load  factors  a  and  b  are  required  to  be 
less  than  2.0  (except  for  tributary  areas  less  than  20  m2 ) . 

3.1.3  Structures  Studied 

An  infinite  number  of  variations  in  types  of  frames 
could  be  studied.  To  limit  the  problem,  symmetrical  two, 
three  and  five  span  frames  based  on  Section  8.9.1  of  the  ACI 
code  (1)  are  studied.  These  structures  are  single-story 
sections  of  frames  with  the  far  ends  of  columns  above  and 
below  the  floor  assumed  to  be  fixed. 

Symmetrical  frames  with  variations  in  the  following 
variables  are  studied: 

a.  the  tributary  area, 

b.  the  number  of  spans, 

c.  the  span  length  ratio, 

d.  the  column  to  beam  stiffness  ratio. 
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Variable  a  affects  the  loads  on  a  span;  variables  b,  c  and  d 
affect  the  influence  coefficient  matrix.  The  frames  studied 
are  listed  in  Table  3.2.  The  values  of  the  variables  are 
chosen  in  order  to  get  the  effects  of  each  variable. 

The  load  effects  of  interest  in  this  study  are  moments 
at  the  ends,  quarter  points  and  mid  points  of  the  spans,  end 
shears,  column  moments  and  axial  forces.  So  as  to  be  able  to 
construct  moment  envelopes,  lifetime  minimum  and  maximum 
values  are  required.  For  the  others,  only  absolute  lifetime 
maximum  values  are  needed.  Since  the  frames  considered 
herein  are  symmetrical,  only  the  load  effects  on  one  half  of 
a  frame  are  required  for  analysis.  Axial  deformations  of 
members  are  ignored  in  the  structural  analysis. 


3.2  Dead  Load  Study 

3.2.1  Statistical  Calculation 

Direct  statistical  calculation  of  dead  load  effects  is 
based  on  the  following  theorem  (3):  If  a  random  variable,  Y, 
is  the  sum  of  normally  distributed  variables,  Xj  ,  namely, 

Y  =  ^  aj  Xj  (3.7) 

i=1 

where  a-  is  the  coefficient  associated  with  X(  ,  then  Y  is 
also  normally  distributed.  The  mean  of  Y  is  (Eq.  2.4.81a  of 
Ref.  3): 

E  (  Y  )  =  a,  E  ( Xj  ) 
i=1 


(3.8) 
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TABLE  3.2  FRAMES  STUDIED 


Frame 

Number  of 
Spans 

Span  Length 
Rati  o 

* 

(m2) 

* 

kb 

1 

2 

1  :1 

25 

0.2 

2 

3 

1  :1  : 1 

25 

0.2 

3 

5 

1  : 1  :1  :1  :1 

25 

0.2 

4 

3 

1  : 1  :1 

25 

2.0 

5 

3 

1  :1  .5:1 

25 

0.2 

6 

3 

1  .5:1  : 1  .5 

25 

0.2 

7 

3 

1  : 1  : 1 

50 

0.2 

8 

3 

1  : 1  :1 

10 

0.2 

* 


Based  on  the  (longer)  span  length  of  5  m 
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and  the  variance  of  Y  is  ( Eq .  2.4.81b  of  Ref.  3): 


Var ( Y) 


af  Var  (  Xj  )  +  2  a j  aj  Cov(  Xj  ,  Xj  ) 


n  0  n-1 


(3.9) 


where  Cov(Xj  , Xj  )  represents  the  covariance  of  variables  X, 
and  Xj  .  The  results  are  valid  whether  X,  are  correlated  or 
independent.  In  the  latter  case  Cov(Xj  , Xj  )  equals  zero  and 
the  second  term  in  Eq.  3.9  disappears. 

A  particular  dead  load  effect,  Yj  ,  is  a  linear  function 
of  the  uniformly  distributed  dead  load  on  each  span.  The 
relation  is  represented  by 


n 


(3.10) 


where  a{-  is  an  element  of  the  influence  coefficient  matrix 
at  row  i  and  column  j  for  a  particular  frame  (derivation  of 
the  matrix  is  shown  in  Appendix  A).  The  term  Xj  is  the 
uniformly  distributed  dead  load  on  span  j,  and  n  is  the 
total  number  of  spans  in  the  frame.  Since  the  dead  loads, 

Xj  ,  are  normally  distributed  variables  (sect.  2.2),  Yj  is 
also  normally  distributed.  The  mean  and  the  variance  of  Xj 
is  the  same  for  j  =  1  to  n,  therefore,  using  Eqs  3.8  and 
3.9,  the  mean  and  variance  of  Yj  become 


n 


E(Y,  )  =  mx  ^3|j 


(3.11) 


(3.12) 
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where  r  is  the  correlation  coefficient  between  Xj  .  As 
mentioned  in  the  dead  load  model  in  Section  2.2,  the  mean 
unit  dead  load  (or  the  mean  uniformly  distributed  load  per 
unit  width),  mx,  is  assumed  to  be  1.05  times  the  nominal 
unit  dead  load  which  is  taken  equal  to  unity  for  analysis, 
and  the  coefficient  of  variation,  Vx ,  is  assumed  to  be  0.10. 
No  data  are  available  to  estimate  the  value  of  the 
correlation  coefficient,  r,  between  Xj  .  However,  it  is  most 
likely  that  the  dead  load  on  one  span  is  highly  correlated 
with  the  dead  load  on  adjacent  spans,  therefore,  r  =  0.75 
has  been  arbitrarily  assumed  for  analysis.  The  extreme 
situations  of  r  =  0.0  (completely  independent  loads  on  all 
spans)  are  also  examined  for  comparison.  Perfectly 
correlated  loads  (r  =  1.0)  have  not  been  considered  since 
this  corresponds  to  constant  loads  on  all  spans.  Based  on 
the  normal  distribution  of  the  dead  load  effect  with  its 
mean  and  variance  given  by  Eqs.  3.11  and  3.12,  any  fractile 
of  a  dead  load  effect  can  be  calculated. 

3.2.2  Linear  Programming 

In  the  similar  manner  as  for  the  live  load,  the  dead 
load  effects  obtained  from  the  statistical  calculation  for 
the  extreme  cases  are  modelled  with  simple  load  patterns. 

The  formulation  of  the  linear  program  for  the  dead  load  is 
similar  to  that  discussed  in  Section  3.1.2(c)  for  the  live 
load.  The  term  Mmj  in  Section  3.1.2(c)  is  redefined  here  as 
the  absolute  values  of  the  0.999  fractiles  determined  from 
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the  statistical  calculation.  For  comparison,  the  load 
pattern  based  on  0.99  fractiles  of  the  dead  load  effects  is 
also  examined. 

3.2.3  Structures  Studied 

The  type  of  structures  studied  for  the  dead  load  is  the 
same  as  for  the  live  load  (Sect.  3.1.3).  Since  the  results 
of  the  analysis  of  live  loads  presented  in  Chapter  4  show 
that  the  type  of  frames  has  an  insignificant  influence  on 
the  final  results  of  analysis,  Frame  2  in  Table  3.2  has  been 
arbitrarily  chosen  to  be  studied  for  the  dead  load. 
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4.  DISCUSSION  OF  RESULTS 


4.1  Live  Load 

4.1.1  Results  of  Monte  Carlo  Analysis 

(a)  Overview  of  Monte  Carlo  Results 

Eight  different  frames  (see  Table  3.2)  have  been 
studied  in  the  Monte  Carlo  analysis.  The  Monte  Carlo  results 
of  a  typical  frame  with  3  spans  of  constant  span  length, 
column  to  beam  stiffness  ratio  equal  to  0.2  and  tributary 
area  of  one  span  equal  to  25  sq .  meters  (Frame  2)  are  shown 
in  Figs.  4. 1-4.3.  Each  figure  shows  the  moment  envelopes 
obtained  from  the  0.999,  0.99  or  0.5  fractiles  for  loading 
Cases  MCI  and  MC3,  and  the  moment  envelopes  from  the  ACI 
load  pattern  (the  traditional  pattern). 

In  Fig.  4.1,  the  0.999  fractiles  of  the  negative 
moments  at  the  interior  support  for  loading  Case  MC3  are 
about  10%  higher  than  those  for  Case  MCI,  while  moment 
envelopes  for  both  cases  are  the  same  at  all  other  sections. 
The  comparison  of  loading  Cases  MCI,  MC2  and  MC3  for 
different  frames  will  be  discussed  in  Section  4.1.1(c). 

Figure  4.1  also  shows  that  the  factored  moment  envelope 
from  the  ACI  load  pattern  based  on  a  load  factor  of  1.7  and 
a  nominal  live  load  given  by  Eq.  3.1  is  close  to  the  0.999 
fractiles  from  the  Monte  Carlo  results,  except  for  negative 
moments  at  the  interior  support  and  positive  moments  at  the 
interior  support.  Since  the  latter  are  offset  by  dead  load 
moments,  they  will  be  disregarded. 
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Fig. 4.1  Comparison  of  0.999  fractiles  from  Monte  Carlo  results 
the  moment  envelope  from  ACI  load  pattern  -  Frame  2. 
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Fig.4.2  Comparison  of  0.99  fractiles  from  Monte  Carlo  results  with 
the  moment  envelope  from  ACI  load  pattern- Frame  2. 
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Fig.4.3  Comparison  of  0.50  fractiles  from  Monte  Carlo  results  with 
the  moment  envelope  from  ACI  load  pattern -Frame  2. 
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The  interior  negative  support  moments  from  the  ACI  load 
pattern  are  about  40%  larger  than  the  0.999  fractiles  for 
Case  MC3.  This  large  discrepancy  suggests  that  the  design 
load  pattern  of  full  factored  load  on  the  two  spans  adjacent 
to  the  support  is  too  severe.  This  is  confirmed  by  the 
results  of  linear  programming  discussed  in  Section  4.1.2. 

In  a  similar  manner  to  Fig.  4.1,  Fig.  4.2  shows  that 
the  0.99  fractiles  of  the  negative  moments  at  the  interior 
support  for  loading  Case  MC3  are  larger  than  those  for 
Case  MCI,  and  the  discrepancy  between  the  0.99  fractiles  of 
the  interior  negative  support  moments  and  those  from  the 
factored  ACI  load  pattern  is  much  larger  than  the 
discrepancy  at  the  other  sections.  In  Fig.  4.3,  the  moment 
envelope  from  0.5  fractiles  is  compared  to  that  from  the 
unfactored  ACI  load  pattern.  The  unfactored  moment  envelope 
from  the  ACI  load  pattern  is  seen  to  be  more  severe  than 
that  from  the  0.5  fractiles  at  all  critical  sections. 

With  respect  to  the  end  shears,  the  0.999  fractiles 
from  the  Monte  Carlo  results  are  slightly  less  than  the 
factored  shears  from  the  ACI  load  pattern.  For  the  interior 
column  axial  load,  the  0.999  fractile  based  on  loading 
Case  MC3  is  about  10%  higher  than  that  based  on  loading 
Case  MCI,  while  the  exterior  column  load  remains  the  same  in 
both  loading  cases.  The  different  loading  cases  will  be 
discussed  more  fully  in  Section  4.1.1(c).  In  the  same  way  as 
the  negative  moment  at  the  interior  support,  the  interior 
column  axial  load  obtained  from  the  ACI  load  pattern  based 
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on  a  load  factor  of  1.7  is  much  larger  than  the  0.999 
fractile.  Again,  this  is  confirmed  by  the  results  of  linear 
programming  discussed  in  Section  4.1.2. 

(b)  Effects  of  Variables  on  Monte  Carlo  Results 

A  number  of  parameters  have  been  varied  in  the  Monte 
Carlo  Analysis.  They  can  be  divided  into  two  categories:  the 
parameters  of  the  vacant  period,  and  those  of  frame 
properties  which  include  the  tributary  area,  number  of 
spans,  span  length  ratio  and  column  to  beam  stiffness  ratio. 

In  the  description  of  the  live  load  model,  it  has  been 
arbitrarily  assumed  that  no  live  loads  are  present  during 
"vacant  periods"  representing  tenant  changes,  etc.  Although 
it  seems  reasonable  that  they  should  occur,  no  data  are 
available  on  the  characteristics  of  such  periods.  In  the 
analysis,  the  vacant  period  has  been  described  in  terms  of 
its  duration  and  the  time  between  successive  vacant  periods. 
To  study  the  effects  of  the  assumed  values,  a  parametric 
study  has  been  made.  Table  4.1  shows  the  effects  of  varying 
the  parameters  of  the  vacant  period  on  the  Monte  Carlo 
results.  The  0.99  fractiles  from  the  500-run  computer 
results  using  Frame  8  (see  Table  3.2)  based  on  the  assumed 
parameters  in  Section  2.1.4  are  compared  with  the  results 
based  on  other  sets  of  parameters.  The  sums  of  the  values, 
while  meaningless  in  themselves,  show  a  small  increase  in 
results  when  the  time  between  periods  increases  or  the 
duration  of  the  period  decreases.  These  differences  are  of 


TABLE  4.1  EFFECTS  OF  VARYING  THE  PARAMETERS  DESCRIBING  THE  VACANT  PERIOD  ON  THE 
MONTE  CARLO  RESULTS  (Based  on  Case  MCI)  -  Frame  8  (See  Table  3.2) 
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Fig. 4.4  Locations  of  Calculated  Load  Effects 
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the  same  magnitude  as  the  errors  introduced  by  other 
assumptions  and  the  values  of  these  parameters  as  chosen  in 
Section  2.1.4  are  assumed  to  be  good  enough  for  the  final 
analysis . 

The  tributary  area  of  a  frame  affects  the  loads  on  each 
span,  while  other  frame  properties  affect  the  influence 
coefficients,  ie,  the  transformation  of  the  loads  to  load 
effects.  In  the  description  of  the  live  load  model,  the 
magnitude  of  the  live  load  is  a  function  of  tributary  area 
(or  influence  area);  the  larger  is  the  tributary  area,  the 
smaller  is  the  live  load.  As  a  result,  for  comparable 
frames,  the  frame  with  a  larger  tributary  area  obtains 
smaller  values  of  load  effects  per  unit  width  from  the  Monte 
Carlo  analysis  than  one  with  a  smaller  tributary  area  does. 
For  the  other  frame  properties,  the  longer  span  of  two 
adjacent  uneven  spans  gets  larger  moments  and  shears  than 
the  shorter  span  does;  when  the  column  to  beam  stiffness 
ratio  is  higher,  the  beam  moments  become  smaller  while  the 
column  moments  are  larger.  The  effects  of  different  number 
of  spans  on  the  Monte  Carlo  results,  however,  are  not 
sufficiently  consistent  to  draw  any  conclusion.  To  account 
for  the  effects  of  a  series  of  frame  variables,  a  number  of 
different  frames  have  been  considered  in  the  linear 
programming  analysis  presented  in  Section  4.1.2(b). 
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(c)  Comparison  of  Monte  Carlo  Results  for  Loading  Cases  MCI, 
MC2  and  MC3 

The  0.999  fractiles  of  the  interior  negative  support 
moments  and  the  interior  column  load  based  on  Cases  MCI  and 
MC2  are  compared  in  Table  4.2  for  four  different  frames.  The 
results  indicate  that  Case  MC2  (two  adjacent  spans  occupied 
by  one  room)  always  govern  the  interior  column  load  but  may 
or  may  not  govern  the  interior  moments.  This  can  be 
explained  by  observing  the  influence  coefficients  in 
Table  4.2.  For  the  interior  column  load,  the  influence 
coefficients  are  always  significant  in  the  first  two 
adjacent  spans  for  any  frame.  For  the  support  moments, 

Case  MCI  governs  when  the  influence  coefficient  in  one  of 
the  first  two  spans  is  relatively  small,  as  shown  by  Frame  4 
in  which  the  column  to  beam  stiffness  ratio  is  2.0. 

The  results  of  the  comparison,  therefore,  verify  the 
assumption  discussed  in  Section  3.1.1(b)  that, 
conservatively,  the  span(s)  where  the  influence  surface  is 
significant  should  be  occupied  by  one  room  in  the  Monte 
Carlo  analysis. 

A  third  loading  case,  Case  MC3,  has  been  derived  using 
the  larger  of  the  Monte  Carlo  values  from  Case  MCI  and 
Case  MC2  for  the  interior  negative  support  moments,  the 
interior  column  loads  from  Case  MC2 ,  and  all  other  moments, 
shears  and  exterior  column  loads  from  Case  MCI.  While  this 
appears  to  be  an  extreme  case,  it  does  allow  for  the 
possibilities  that  the  loads  on  any  two  adjacent  spans  may 
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TABLE  4.2  COMPARISON  OF  MONTE  CARLO  RESULTS  BASED  ON  CASES  MCI  AND  MC2 
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be  correlated  or  that  the  loads  on  all  spans  are 
i ndependent . 

4.1.2  Results  of  Linear  Programming 

(a)  Comparison  of  Different  Design  Load  Patterns 
The  design  load  patterns  considered  in  this  study  are 
based  on  three  types  of  design  loading  illustrated  in 
Fig.  4.5.  These  have  been  selected  because  of  their 
simplicity  and  their  ability  to  represent  extreme  loading 
cases.  The  three  types  of  loading  are: 

alternate  span  loading  -  heavy  factored  design  load, 

a1  Ln ,  on  alternate  spans,  and  light  factored  design 
load,  b1 Ln  ,  on  all  other  spans  where  a1  >  b1  , 

single  span  loading  -  heavy  factored  design  load,  a2l_n, 
on  one  span  and  light  factored  design  load,  b2Ln,  on 
all  other  spans  where  a2  >  b2 , 

adjacent  span  loading  -  heavy  factored  design  load, 
a3 Ln ,  on  two  adjacent  spans  and  light  factored 
design  load,  b3Ln ,  on  all  other  spans  where  a3  >  b3 . 

In  each  type  of  loading,  the  load  on  a  span  is  the  product 
of  a  load  factor,  a  or  b,  and  the  nominal  live  load,  Ln , 
which  is  given  by  Eq .  3.1. 

From  the  three  basic  types  of  loading,  four  design  load 
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Fig.  4.5  Types  of  Design  Loading 
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patterns  are  considered: 

Pattern  A  -  a  combination  of  alternate  span  loading 

(Fig.  4.5a)  plus  adjacent  span  loading  (Fig.  4.5c) 
as  appropriate,  with  a1 =  a3  and  b1 =  b3 , 

Pattern  B  -  single  span  loading  (Fig.  4.5b) 

Pattern  C  -  a  combination  of  alternate  span  loading 

(Fig,  4.5a)  plus  adjacent  span  loading  (Fig.  4.5c) 
as  appropriate,  with  a1  and  b1  independent  of  a3  and 
b3  respectively, 

Pattern  D  -  a  combination  of  single  span  loading 

(Fig.  4.5b)  plus  adjacent  span  loading  (Fig.  4.5c) 
as  appropriate,  with  a2  and  b2  independent  of  a3  and 
b3  respectively. 

Pattern  A  is  similar  to  the  traditional  design  load  pattern 
except  that  there  may  be  a  small  load  on  the  unloaded  spans. 
Pattern  B  is  a  simpler  model  suggested  by  Beeby  (2).  The 
consideration  of  Patterns  C  and  D  results  from  the  lack  of 
consistency  in  the  results  of  Patterns  A  and  B  as  discussed 
later.  A  possible  load  pattern  similar  to  Pattern  D  except 
a2  =  a3  and  b2  =  b3  is  not  considered  because  the  results  of 
Pattern  D  indicates  that  a2  is  much  different  from  a3  as 


shown  later. 
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The  0.999  fractiles  of  the  lifetime  minimum  and  maximum 
load  effects  from  the  Monte  Carlo  analysis  of  Frame  2  (with 
constant  span  length  ratio  and  column  to  beam  stiffness 
ratio  equal  to  0.2,  see  Table  3.2)  have  been  used  to 
investigate  the  validity  of  different  design  load  patterns. 
The  results  are  presented  in  Table  4.3. 

Pattern  A 

Based  on  loading  Case  MCI  (loads  on  all  spans 
independent),  load  factors  of  1.56  and  0.02  are  obtained  for 
heavily  and  lightly  loaded  spans.  The  penalty  obtained  is 
very  high  because  the  design  load  pattern  overestimates  the 
interior  negative  support  moments  by  more  than  40%.  In  the 
design  load  pattern,  those  load  effects  are  modelled  with 
the  adjacent  span  loading.  The  other  load  effects  which  are 
modelled  with  the  alternate  span  loading  are  relatively  well 
f i t  ted . 

Loading  Case  MC3  gives  higher  values  for  the  interior 
support  moments  and  the  interior  column  load  than  Case  MCI 
does  (see  Table  4.2).  When  these  load  effects  are  modelled 
with  Pattern  A,  better  agreement  is  obtained.  Nevertheless, 
they  are  still  significantly  overestimated,  as  shown  by  the 
high  penalty.  These  results,  however,  infer  that  the  penalty 
can  be  minimized  by  treating  load  factors  of  the  alternate 
span  loading  and  those  of  the  adjacent  span  loading 
independently,  leading  to  consideration  of  Pattern  C. 
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TABLE  4.3  COMPARISON  OF  DIFFERENT  DESIGN  LOAD  PATTERNS 


Pattern  A 

Pattern  B 

Pattern  C 

Pattern  D 

Load 

Factor 

Case 

MCI 

Case 

MC3 

Case 

MCI 

Case 

MC3 

Case 

MCI 

Case 

MC3 

Case 

MCI 

Case 

MC3 

al  ,2 

1.56 

1.56 

1  .66 

1  .66 

1  .63 

1  .63 

1  .66 

1  .66 

b1  >2 

0.02 

0.02 

0.07 

0.07 

0.06 

0.06 

0  .03 

0  .03 

a3 

- 

- 

- 

- 

1.26 

1  .24 

1  .26 

1  .24 

b3 

- 

- 

- 

- 

1  .23 

0.28 

1  .23 

0.28 

Total 
Penal ty 

12.6 

10.5 

7.3 

11.4 

4.7 

4.0 

2.9 

2.2 
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Pattern  B 

When  load  factors  are  derived  for  Pattern  B  based  on 
loading  Case  MCI,  the  penalty  obtained  is  less  than  for 
Pattern  A.  For  Case  MCI,  Pattern  B  underestimates  the 
interior  support  moments  and  the  interior  column  load  by 
more  than  10%,  while  other  load  effects  are  well  fitted.  For 
Case  MC3,  the  penalty  is  even  higher  because  the  interior 
support  moments  and  the  interor  column  load  are  more 
severely  underestimated.  This  is  conceivable  because  for 
those  load  effects  with  significant  influence  surface  in 
adjacent  spans,  the  single  span  loading  pattern  (Pattern  B) 
appears  less  representative  of  the  extreme  loading  cases 
than  the  adjacent  span  loading  does.  Therefore,  it  is 
reasonable  to  also  consider  a  Pattern  D,  which  includes  the 
adjacent  span  loading  in  addition  to  the  single  span 
loading . 

Patterns  C  and  D 

The  load  factors  for  the  Pattern  C  loading  are  given  in 
Table  4.3.  These  suggest  that  the  0.999  fractiles  of  the 
load  effects  from  the  Monte  Carlo  analysis  based  on  Case  MCI 
can  best  be  represented  by  an  alternate  span  loading  with 
1.63  or  0.06  times  the  nominal  load  placed  to  give  maximum 
values  of  load  effects  normally  predicted  by  alternate  span 
values,  plus  a  constant  load  of  1.23  to  1.26  times  the 
nominal  load  when  predicting  interior  column  loads  or 
negative  moments  at  interior  supports.  For  Case  MC3,  the 
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load  factors  of  the  alternate  span  loading  are  the  same  as 
those  of  Case  MCI  because  the  alternate  span  loading  fits 
the  same  load  effects  in  both  loading  cases.  The  adjacent 
span  loading,  however,  is  different,  with  1.24  times  the 
nominal  load  on  two  adjacent  spans  and  a  small  load  on  the 
other  spans.  The  penalty  for  Pattern  C,  shown  in  Table  4.3, 
is  considerably  reduced.  As  anticipated,  the  load  factor  a3 
of  the  adjacent  span  loading,  which  models  the  interior 
support  moments  and  the  interior  column  load,  is 
significantly  less  than  the  load  factor  a1  of  the  alternate 
span  loading. 

In  Pattern  D,  the  load  factors  for  the  adjacent  span 
loading  are  the  same  as  those  in  Pattern  C  because  they  both 
model  the  same  load  effects,  while  the  load  factors  for  the 
single  span  loading  are  slightly  higher  than  those  in 
Pattern  D.  The  penalty  for  Pattern  D  is  less  than  for 
Pattern  C,  which  indicates  that  the  single  span  loading  may 
be  more  representat i ve  of  the  extreme  loading  cases  for 
0.999  fractiles  of  the  load  effects  concerned  than  the 
alternate  span  loading. 

Patterns  C  and  D  wi 1 1  be  examined  in  the  next  section 
for  different  frames  to  investigate  the  effects  of  other 
variables  and  also  as  a  continuous  comparison  of  the  two 
load  patterns.  Because  the  adjacent  span  loading  models  the 
same  load  effects  in  both  patterns,  the  real  comparison  is 
between  the  alternate  span  and  single  span  loadings,  which 
model  the  load  effects  other  than  those  modelled  by  the 
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adjacent  span  loading.  The  significance  of  the  penalty  as 
related  to  the  degree  of  fit  of  a  pattern  will  be 
illustrated  as  well  in  the  next  section. 

(b)  Effects  of  Variables 

The  load  factors  of  Patterns  C  and  D  for  different 
frames  based  on  Cases  MCI  and  MC3  are  presented  in 
Table  4.4.  As  mentioned  previously,  the  load  factors  of  the 
adjacent  span  loading  are  the  same  in  both  Patterns  C  and  D. 
With  two  exceptions  discussed  later,  the  adjacent  span 
loading  has  been  used  to  fit  the  interior  negative  support 
moments  and  the  interior  column  loads.  The  other  load 
effects  have  been  modelled  with  the  single  span  or  alternate 
span  loading.  The  load  factors  for  the  single  span  or 
alternate  span  loading  are  essentially  the  same  in  Cases  MCI 
and  MC3,  because  in  both  cases  the  single  span  or  alternate 
span  loading  has  been  used  to  fit  the  same  load  effects. 

It  has  been  shown  in  Table  4.2  that  Case  MCI  governs 
the  interior  negative  support  moments  in  Frame  4  (with 
column  to  beam  stiffness  ratio  equal  to  2.0,  see  Table  3.2) 
and  also  governs  the  interior  negative  support  moment  of  the 
longer  (exterior)  span  in  Frame  6  (with  span  length  ratio 
equal  to  1.5:1. 0:1. 5,  see  Table  3.2).  This  is  because  the 
influence  surface  in  one  of  the  adjacent  spans  is  relatively 
small,  as  explained  in  Section  4.1.1(c).  When  considering 
Frame  4,  the  interior  support  moments  and  the  interior 
column  load  have  been  poorly  fitted  by  the  adjacent  span 
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TABLE  4.4  EFFECTS  OF  VARIABLES  ON  DESIGN  LOAD  PATTERNS 
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loading;  the  penalty  is  severe  because  the  interior  column 
load  is  much  overestimated  while  the  negative  moments  are 
underestimated  if  a  constant  load  factor  a3  is  used.  The 
load  factor  a3  is  about  1.35,  which  is  higher  than  the  value 
obtained  by  Frame  2  presented  in  Table  4.3.  To  reduce  the 
penalty  and  hence  increase  the  accuracy,  the  interior 
negative  support  moments  have  been  included  in  the  single 
span  or  alternate  span  loading,  which  gives  much  better  fit, 
especially  the  single  span  loading.  The  adjacent  span 
loading,  which  fits  the  interior  column  load  only,  gives  the 
load  factor  a3  equal  to  1.07  for  Case  MCI  and  1.20  for 
Case  MC3,  as  presented  in  Table  4.4. 

In  a  similar  manner,  the  interior  support  moment  of  the 
longer  span  in  Frame  6  is  also  included  in  the  single  span 
or  alternate  span  loading.  The  above  special  arrangement  of 
the  interior  support  moments,  however,  leads  to  the 
consideration  of  both  single  span  (or  alternate  span) 
loading  and  adjacent  span  loading  even  when  calculating  the 
interior  negative  support  moments. 

Table  4.4  shows  that  the  load  factors  of  the  alternate 
span  loading  or  the  single  span  loading  appear  insensitive 
to  most  variables  except  the  tributary  area  (Frames  7,  8). 
For  a  tributary  area,  AT ,  of  10  m2(Frame  8),  the  load  factor 
a1  or  a2  is  much  higher  than  the  load  factor  for  AT  =  25  m 
(Frames  1-6).  This  is  because  the  nominal  live  load 
(Eq.  3.1)  becomes  constant  for  AT  <  20  m2 ,  while  the  loads 
determined  from  the  live  load  model  as  used  in  this  study 
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Keeps  increasing  for  AT  <  20  m2(see  Eq.  2.7,  Fig.  2.1).  For 
At  =  50  m2 ,  the  load  factor  a  is  smaller  than  that  for  AT 
=  25  m2 ,  but  the  effect  can  be  considered  relatively 
insignificant.  In  all  cases,  the  load  factor  b1  or  b2  is 
close  to  zero. 

The  variables  have  some  significant  effects  on  the 
penalties  in  the  alternate  span  loading  as  shown  in 
Table  4.4.  The  penalties  in  the  single  span  loading, 
however,  consistenly  remain  small  and  less  than  the 
penalties  in  the  alternate  span  loading  for  all  different 
frames.  Note  that  the  penalty  for  the  five  span  frame 
(Frame  3)  is  higher  than  the  others  because  more  load 
effects  have  been  fitted.  This  indicates  that  the  single 
span  loading  is  consistently  a  better  representation  of  the 
extreme  loading  cases  for  the  load  effects  with  a 
significant  influence  surface  in  one  span  and  relatively 
small  influence  surfaces  in  all  other  spans.  The  average 
value  of  the  load  factor  a2  of  the  single  span  loading  for 
Frames  1-7  is  1.60  (excluding  the  large  load  factor  obtained 
by  Frame  8),  and  b2  is  practically  equal  to  zero. 

In  a  similar  manner  as  for  the  alternate  span  or  single 
span  loadings,  most  variables  have  an  insignificant  effect 
on  the  load  factors  for  the  adjacent  span  loading  case 
except  the  effect  of  small  tributary  area.  The  load  factor 
a3  is  similar  for  Cases  MCI  and  MC3  ranging  from  1.07  to 
1.26  for  MCI  and  1.20  to  1.26  for  MC3  for  frames  with  large 
tributary  areas  (Frames  1-7).  Since  Case  MC3  represents  a 
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more  conservative  situation,  it  is  considered  more  desirable 
as  the  basis  for  a  design  recommendation.  The  values  of  the 
load  factor  a3  for  different  frames  (except  Frame  8)  based 
on  Case  MC3  are  consistent  and  the  average  is  1.23,  which  is 
0.77  of  the  load  factor  a2  of  the  single  span  loading.  For 
simplicity,  it  is  proposed  that  in  design  the  load  factor  a3 
should  be  taken  equal  to  3/4  of  the  load  factor  of  the 
single  span  loading.  The  load  factor  b3 ,  which  ranges  from  0 
to  0.28  for  Frames  1  to  7  with  an  average  of  0.07,  shall  be 
set  equal  to  zero  for  design  purposes. 

The  degree  of  fit  of  Pattern  D  is  illustrated  in 
Tables  4.5  and  4.6  by  Frame  2  (Table  3.2)  for  loading 
Case  MC3  with  load  factors  presented  in  Table  4.4.  Column  5 
in  these  tables  illustrates  the  relative  significance  of  the 
penalty  as  functions  of  a  single  beam  moment,  shear,  etc. 
Pattern  D,  based  on  Case  MC3,  is  taken  as  a  design 
recommendation  which  will  be  presented  more  fully  in  the 
next  section. 

4.1.3  Des i gn  Recommenda  t i on 

It  is  recommended  that  the  design  value  of  a  live  load 
effect  in  a  continuous  beam  or  a  one  way  slab  be  determined 
using  whichever  of  the  following  types  of  loading  that 
produces  the  largest  value: 

a.  factored  live  load  on  one  span, 

b.  3/4  of  the  factored  live  load  on  two  adjacent  spans. 
The  load  factor  correspondi ng  to  0.999  fractiles  of  the 
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TABLE  4.5  COMPARISON  OF  0.999  FRACTILES  OF  LIFETIME  MAXIMUM 
LIVE  LOAD  EFFECTS  WITH  THOSE  FROM  LOAD  PATTERN  D 
-  FRAME  2  (See  Table  3.2) 


Locati on 
(Fig. 4. 4) 

[1] 

0.999 

Fracti 1 es 

[2] 

Desi gn 

Val  ue 

[3] 

Devi  ati  on 
=Abs [ 3]- Abs [2] 

[4] 

■Lil  x  1  00% 

[5] 

1 

0.505 

0.575 

0.070 

1  .0 

2 

4.864 

4.863 

-0.001 

-0.01 

3 

6.847 

6.859 

0.012 

0.2 

4 

5 

3.100 

3.100 

0.0 

0.0 

6 

7 

2.862 

3.127 

0.265 

3.8 

8 

5.642 

6.004 

0.362 

5.2 

9 

8.503 

8.503 

0.0 

0  .0 

10 

- 

- 

- 

- 

11 

9  .078 

9.208 

0.130 

2.3 

12 

1  .442 

1.444 

0.002 

0.04 

13 

0.977 

0.974 

-0.003 

-0.06 

14 

8.502 

8.504 

0.002 

0.04 

★ 

1  5 

15.518 

15.536 

0.018 

0.16 

*  Fitted  by  the  adjacent  span  loading  (see  Sect.  4.1.2(a)  ) 


**  For  beam  moments,  F  = 


LnS2 

—5—  where  L  is  the  nominal  live 
o  n 


load  given  by  Eq.3.1  and  S  is  the 
span  length  (5m). 

LnS 


For  end  shears,  F  =  — 9 

LnS 

For  column  moments,  F  =  — sy 

L  S 

For  column  loads,  F  =  —5—  (exterior) 


F  =  LnS  (interior) 


I  }'  o 

» 
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TABLE  4.6  COMPARISON  OF  0.999  FRACTILES  OF  LIFETIME  MINIMUM 
LIVE  LOAD  EFFECTS  WITH  THOSE  FROM  LOAD  PATTERN  D 
-  FRAME  2  (See  Table  3.2) 


* 

* 

* 

* 

* 

[1] 

[2] 

[3] 

[4] 

[5] 

1 

-2.883 

-2.888 

0.005 

0.07 

2 

-0.347 

-0.337 

-0.010 

-0.1 

3 

-1.200 

-1 .365 

0.165 

2.4 

4 

-2.255 

-2.510 

0.255 

3.7 

5** 

-7.326 

-7.326 

0.0 

0  .0 

** 

6 

-7.107 

-7.107 

0.0 

0.0 

7 

-2.771 

-2.884 

0.113 

1  .6 

8 

-1.787 

-1 .420 

-0.367 

-5.3 

9 

- 

- 

- 

- 

10 

-9.941 

-9.913 

-0.028 

-0.5 

11 

- 

- 

- 

- 

*  Re  f  e  r 

**  Fitted 

to  Table  4.5 

by  the  adjacent 

span  loading 

(see  Sect. 

4.1  .2(a)  ) 

■0 
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lifetime  minimum  and  maximum  live  load  effects  for  office 
buildings  is  1.6.  This  design  load  pattern  is  valid  for 
calculating  moment  envelopes,  end  shears,  column  moments  and 
column  loads. 

The  one  span  loading  gives  larger  values  for  most  load 
effects  than  the  adjacent  span  loading  does.  The  adjacent 
span  loading,  however,  gives  larger  values  for  load  effects 
with  significant  influence  surfaces  in  two  adjacent  spans 
such  as  the  interior  column  loads  and  the  negative  moments 
at  interior  supports.  Nevertheless,  for  the  frame  with  a 
high  ratio  of  column  to  beam  stiffness  ratio  such  as  2.0, 
the  one  span  loading  also  governs  the  negative  moments  at 
interior  supports. 


4.2  Dead  Load 

The  British  Standard  Code  of  Practice  (4)  requires 
pattern  application  of  dead  load  to  a  structure  as  shown  in 
Table  4.7.  This  is  not  required  by  the  ACI  Code  (1)  or  NBC 
Code  (19).  The  0.999  and  0.99  fractiles  of  the  dead  load 
effects  have  been  found  to  be  well  modelled  by  Load 
Pattern  A  (see  Sect.  4.1.2(a)),  which  is  similar  to  the 
traditional  load  pattern.  The  load  factors,  a  and  b, 
corresponding  to  the  0.999  and  0.99  fractiles  are  presented 
in  Tables  4.8  and  4.9  for  three  values  of  the  coefficient  of 
correlation,  r,  between  the  dead  loads  on  adjacent  spans. 

Table  4.8  shows  that  the  British  dead  load  pattern 
closely  approaches  the  load  pattern  which  produces  0.999 
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TABLE  4.7  CURRENT  DESIGN  DEAD  LOAD  FACTORS  STIPULATED 
BY  DIFFERENT  DESIGN  CODES 


Code 

Bri  tish 

NBC 

ACI 

* 

a 

1  .4 

1  .25 

1  .4 

** 

b 

1  .0 

1  .25 

1  .4 

TABLE  4.8  LOAD  FACTORS  CORRESPONDING  TO  0.999  FRACTILES 
OF  DEAD  LOAD  EFFECTS 


Correl ati on 
Coeffi  ci  ent 

0.0 

0.75 

1  .0 

* 

a 

1  .35 

1.37 

1.37 

** 

b 

0.96 

1  .22 

1.37 

Total 

Penal ty 

0.86 

0.35 

- 

TABLE  4.9  LOAD  FACTORS  CORRESPONDING  TO  0.99  FRACTILES 
OF  DEAD  LOAD  EFFECTS 


Correl ati on 
Coeffi  ci  ent 

0.0 

0.75 

1  .0 

* 

a 

1  .27 

1  .29 

1  .29 

** 

b 

0.98 

1  .18 

1  .29 

Total 

Penalty 

0.65 

0.26 

- 

*  Large  load  factor  in  Pattern  A  (see  Sect.  4.1.2(a)  ) 
**  Small  load  factor  in  Pattern  A  (see  Sect.  4.1.2(a)  ) 
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fractiles  of  load  effects  for  r  =  0.0  or  completely 
uncorrelated  dead  loads  on  adjacent  spans.  On  the  other 
hand,  the  NBC  and  ACI  assumptions  of  constant  factored  dead 
loads  on  all  spans  approach  the  case  of  strongly  or 
completely  correlated  dead  loads  on  all  spans.  The  dead  load 
factors  in  the  NBC  and  ACI  Codes  correspond  to  the  0.99  and 
0.999  fractiles,  respectively,  as  shown  in  Tables  4.8  and 
4.9. 

The  extreme  situation  of  r  =  0.0  is  unlikely  in  a  real 
situation  since  the  same  source  of  concrete  will  generally 
be  used  for  all  spans  and  the  same  framing  crew  will 
probably  build  the  forms  for  all  spans.  A  more  realistic 
situation  of  r  =  0.75,  however,  indicates  in  Tables  4.8  and 
4.9  that  the  dead  loads  are  approximately  constant  on  all 
spans  for  0.99  and  even  for  0.999  fractiles  of  load  effects. 
For  computational  simplicity,  it  is  reasonable  to  consider 
that  the  factored  dead  loads  are  constant  on  all  spans. 


5.  SUMMARY  AND  CONCLUSIONS 


5.1  Live  Load 

A  Monte  Carlo  process  was  used  to  generate  extreme 
values  of  moments  and  other  major  load  effects  for  a  series 
of  frames  subjected  to  varying  live  loads  during  their 
lifetime.  The  live  load  model  used  in  this  study  was 
developed  by  several  investigators  (11,16,20)  based  on  the 
office  live  loads.  The  analyses  generated  load  effects  for 
loadings  assumed  to  be  independent  on  each  span  and  for 
loadings  which  were  correlated  on  two  adjacent  spans  but 
independent  on  all  others. 

Linear  programming  was  used  to  investigate  the  validity 
of  several  possible  design  load  patterns  for  approximatig 
the  0.999  fractiles  of  lifetime  minimum  and  maximum  load 
effects  from  the  Monte  Carlo  results.  The  results  indicate 
that  the  design  load  pattern  which  best  represents  the  0.999 
fractiles  is  either  the  factored  live  load  on  one  span,  or 
3/4  of  the  factored  live  load  on  any  two  adjacent  spans,  as 
appropriate  to  give  maximum  load  effects.  The  corresponding 
load  factor  was  found  to  be  1.6  to  represent  the  0.999 
fractile  of  the  maximum  loading  effects.  This  design  load 
pattern  represents  a  conservative  loading  assumption  that 
the  loads  may  be  independent  on  all  spans  or  correlated  on 
any  two  adjacent  spans. 

This  design  load  pattern  was  shown  to  be  valid  for 
different  frame  properties  although  when  the  tributary  area 
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of  one  span  is  less  than  20  square  meters,  the  load  factors 
should  be  increased.  This  is  because  the  nominal  live  load 
( Eq .  3.1)  is  unconservative  for  tributary  areas  less  than  20 
square  meters. 


5 . 2  Dead  Load 

Direct  statistical  calculation  of  dead  load  effects  for 
the  extreme  cases  was  performed  due  to  the  simplicity  of  the 
dead  load  model.  Again,  linear  programming  was  used  to 
develop  an  equivalent  design  dead  dead  load  pattern.  The 
results  indicate  that  the  effect  of  pattern  loadings  is 
insignificant  for  a  correlation  coefficient  between  dead 
loads  on  adjacent  spans  equal  to  0.75,  Hence  it  is 
recommended  that  pattern  dead  loads  not  be  used  in  design. 


71 


REFERENCES 


1.  ACI  Standard  318-77,  "Building  Code  Requirements  for 

Reinforced  Concrete  (ACI  318-77)",  American  Concrete 
Institute,  1977. 


2.  Beeby,  A.W.,  "Analysis  of  Slabs",  Cement  and  Concrete 
Association,  Slough,  England  (In  Press). 


3.  Benjamin,  J.R.  and  Cornell,  C.A.,  "Probability 

Statistics  and  Decision  for  Civil  Engineers",  McGraw 
Hill  Book  Co.,  New  York,  1970. 


4.  British  Standard  Code  of  Practice  for  the  Structural  Use 
of  Concrete,  CP1 10,  British  Standards  Institution, 
1972. 


5.  Chalk,  P.L.,  "A  Probabilistic  Determination  of  Design 
Live  Loads  for  Buildings",  Thesis  Submitted  to 
Northwestern  University  in  Partial  Fulfilment  of  the 
Requirements  for  the  Degree  of  Master  of  Science, 
June  1979. 


6.  Corot  is,  R.B.  and  Doshi ,  V.A.,  "Probability  Models  for 

Live  Load  Survey  Results",  Journal  of  the  Structural 
Division,  Vol .  103,  No.  ST6,  June  1979, 
pp.  1257-1274. 


7.  Corotis,  R.B.  and  Jarvia,  V.,  "Stochastic  Considerations 
for  Office  Floor  Loads",  ASCE  Specialty  Conference 
on  Probabilistic  Mechanics  and  Structural 
Probability,  Tucson,  Ar . ,  Januaury  1979. 


8.  Cross,  H.  and  Morgan,  N.D.,  "Continuous  Frames  of 

Reinforced  Concrete",  Chapter  6,  John  Wiley  &  Sons 
Inc.,  1932. 


9.  Culver,  C.G.,  "Survey  Results  for  Fire  Loads  and  Live 
Loads  in  Office  Buildings",  NBS  Building  Science 
Series  85,  National  Bureau  of  Standards,  Washington, 
D.C. ,  May  1976. 


72 


10.  Daellenbach,  H.G.  and  Bell,  E.J.,  "User's  Guide  to 
Linear  Programming",  Prentice  Hall  Inc., 
Englecl i f fs ,  N  .  J  .  ,  1970 . 


11.  Ellingwood,  B.R.  and  Culver,  C.G.,  "Analysis  of  Live 
Loads  in  Office  Buildings",  Journal  of  the 
Structural  Division,  ASCE,  Vol .  103,  No.  ST8,  August 
1977. 


12.  Ellingwood,  B.R.,  Galambos,  T.V.,  MacGregor,  J.G., 

Cornell,  C.A.,  "Development  of  a  Probability  Based 
Load  Criterion  for  American  National  Standard  A58 
Building  Code  Requirements  for  Minimum  Design  Loads 
in  Buildings  and  Other  Structures" ,  NBS  Special 
Publication  577,  National  Bureau  of  Standards,  April 
1980. 


13.  Hauser,  R.,  "Load  Correlation  Models  in  Structural 

Reliability",  Thesis  Presented  to  the  Massachusettes 
Institute  of  Technology  in  Partial  Fulfilment  of  the 
Requirements  for  The  Degree  of  Master  of  Science, 
1971. 


14.  IMSL  Library  Reference  Manual,  Edition  7,  IMSL  Inc., 
Texas,  Januaury  1979. 


15.  Lind,  N.  and  Nowak,  A.,  "Risk  Analysis  Procedure  (I)", 
Department  of  Civil  Engineering,  University  of 
Waterloo,  Waterloo,  Ontario,  Canada,  1977. 


16.  McGuire,  R.K.  and  Cornell,  C.A.,  "Live  Load  Effects  in 
Office  Buildings",  Research  Report  R 7 3 - 2 8 , 
Massachusetts  Institute  of  Technology,  Department  of 
Civil  Engineering,  Cambridge,  Mass.,  May  1973. 


17.  McGuire,  R.K.  and  Cornell,  C.A.,  "Live  Load  Effects  in 
Office  Buildings",  Journal  of  the  Structural 
Division,  ASCE,  Vol.  100,  No.  ST7,  July  1974, 
pp.  1351-1366. 


18.  Mitchell,  G.R.  and  Woodgate,  R.W.,  "Floor  Loadings  in 

Office  Buildings  -  The  Results  of  a  Survey",  CP3/71, 
Building  Research  Station,  United  Kingtom,  Januaury 
1971. 


73 


19.  National  Building  Code  of  Canada  (NBC),  by  the  Associate 
Committee  on  the  National  Building  Code,  National 
Research  Council,  Ottawa,  Canada,  1980. 


20.  Peir,  J.C.,  "A  Stochastic  Live  Load  Model  for 

Buildings",  Research  Report  R 7 1  - 3 5 ,  Massachusetts 
Institute  of  Technology,  Department  of  Civil 
Engineering,  Cambridge,  Mass.,  September  1971. 


21.  Peir,  J.C.  and  Cornell,  C.A.,  "Spatial  and  Temporal 
Variability  of  Live  Loads",  Journal  of  the 
Structural  Division,  ASCE,  Vol .  99,  No.  ST5,  May 
1973,  pp.  903-922 


74 


APPENDIX  A 

Formulation  of  an  Influence  Coefficient  Matrix 


Load  effects  caused  by  any  uniformly  distributed  load 
on  each  span  are  represented  by  the  following  relation: 


n 


where 


Y: 


m 

n 


value  of  load  effect  i 

an  element  of  the  influence  coefficient  matrix 
at  row  i  and  column  j  for  a  particular  frame 
uniformly  distributed  load  on  span  j 
total  number  of  load  effects  considered 
total  number  of  spans  in  the  frame 


Accordingly,  are  determined  by  the  following  method: 

Let  XJ=k  =  1  and  X^k  =  0  (A. 2) 

Eq.  A. 2  is  then  applied  in  Eq .  A.1,  which  results  in 

a  ,  =  Y.  for  i  =  1  ,m  ( A .  3 ) 

I  K  I 

because  Yj  =  1  m  can  be  determined  (when  XJ  =  k  =  1  and  X^  k  =0) 
using  a  simple  frame  analysis  computer  program,  the 
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influence  coefficient  matrix  can  be  formulated  by  repeating 
the  analysis  for  K  =  1  to  n. 


APPENDIX  B 


Listing  of  the  Monte  Carlo  Program 


1 

2 

3 

4 

5 

6 

7 

6 

9 

1  O 

1  1 

1  2 

1  3 

1  4 

1  5 

1  6 

1  7 

1  S 

1  9 

20 

2  1 

22 

23 

24 

25 

2  6 

27 

26 

29 

30 

3  1 

22 

33 

34 

35 

3  S 

37 

38 

39 

4  0 

4  1 

42 

43 

46 

45 

4  6 

47 

4  8 

4  9 

50 

S  1 

52 

53 

54 

55 

5  6 

57 

56 

59 

60 

6  1 

62 

63 

66 

65 

66 

67 

66 

69 

70 

7  1 

72 

73 

74 

75 

7  6 

77 

78 

79 

80 

8  1 

82 

83 

84 

85 

86 

8  7 

88 

89 

90 

9  1 

92 

93 

94 

95 

9  6 

9  7 

9  6 

9  9 

1  00 

1  0  1 

102 

103 

1  04 

1  05 

106 

107 

1  06 
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r****»*s****s*s«»**a 


MONTEL  I  VE 


l**»*S**B»S»*»*»*SS*S**SSSS»* 


c** 

e»* 

c**  * 

c**  * 

C**  THIS  PROGRAM  STATISTICALLY  COMPUTES  THE  LIFETIME  MAXIMUM  AND  * 

C**  MINIMUM  LIVE  LOAD  (GRAVITY  LOAD)  EFFECTS  ( MO ME N T S . S H E A RS  AND  • 

C**  AXIAL  FORCES)  IN  A  CONTINUOUS  FRAME  USING  THE  MONTE-CARLO  * 

C**  TECHNIOUE  • 

C*  *  * 

c 

BAS  I C  ASSUMPT I ONS  : 

LIVE  LOADS  ARE  COMPOSED  OF  SUSTAINED  LOADS  AND  TRANSIENT 
LOADS.  PERIODS  OF  UNOCCUPANCY  ARE  ALSO  CONSIDERED.  THE 
LOADING  ON  ONE  SPAN  IS  UNI  FORM  (OVER  AN  AREA  EOUAL  TO 
TWICE  THE  TRIBUTARY  AREA  OF  ONE  SPAN)  AND  STATISTICALLY 
INDEPENDENT  OF  THE  LOADING  ON  ANY  OTHER  SPAN.  SUCCESSIVE 
LOADING  EVENTS  ARE  ALSO  INDEPENDENT. 


REMARKS : 

THE  NUMBER  OF  RUNS  (NURUN)  MUST  BE  GREATER  THAN  1. 

THE  NUMBER  OF  SPANS  MUST  NOT  BE  LESS  THAN  2. 

IF  VALUES  OF  0.01  FRACTILE  TO  0.99  FRACTILE  WITH  INCREMENTS 
OF  0.01  FRACTILE  ARE  DESIRED  IN  THE  OUTPUT, NURUN  MUST  BE 
MULTIPLES  OF  100 

DOUBLE  PRECISION  S S E E D , S D S E E D , U NS E E D , T L S E E D . T S E E D 
COMMON  SSEED , SDSEED ,  UNSEEO  ,  TLSEED , TSEED 

COMMON  XMAXMl  5000, 23  )  ,XMINM(5000,  18)  ,0(500,5)  ,  INT(S00,5)  . 

1  TRAN  L  (  200 )  ,  1 NTR  | 200 )  , S  US  L  N (  100)  ,  1  NTSN  (  100)  ,XMC0EF(23,5)  , 

1  PS ( 5  )  ,  P  T  ( 5  )  ,  S  L  ENDA  (  5  )  ,TLENDA(5)  ,  1SL I FE, NURUN, NUSPAN, NUROW, 

1  ISMEAN, ISDMIN,NDMIN,NDMAX,NTIME,NFMIN, ITMEAN, ITRMIN 

COMMON  ISPAN,  I  RUN  ,  NUO 1  .NFRAME,  NUR0W2 

CALL  DATA 

DO  100  I  RUN  :  1  ,  NURUN 
DO  200  I  SPAN  *  1  , NUSPAN 
CALL  SUSTAN 
CALL  TRANS'! 

CALL  COMB  N  D 
CONTINUE 
CALL  EFFECT 


1  00 
C 

c 

c 


CONTINUE 
CALL  OUTPUT 
STOP 


END 

C 

c 

c  ***************** 

c 

C  SUBROUTINE  DATA 

C 

C  SSSX*SS*SX*«***** 

c 

C 

C  THIS  SUBROUTINE  READS  AND  WRITES  ALL  THE  INPUT  DATA  AND  THEN 
C  CALCULATES  THE  REQUIRED  STATISTICAL  PARAMETERS.  ALL  THE  TIME 
C  RELATED  VALUES  ARE  CONVERTED  INTO  UNITS  OF  TRANSIENT  LOAD 

C  DURAT I  ON  , THAT  IS.ONE  TIME  UNIT  IS  EOUAL  TO  THE  TIME  OF  THE 

C  DURATION  OF  ONE  TRANSIENT  LOAD.  ALL  SUBSEQUENT  CALCULATIONS 
C  WILL  BE  BASED  ON  THIS  UNIT. 

C 

SUBROUTINE  DATA 


C 


DOUBLE 


COMMON 

COMMON 


1 

1 

1 


COMMON 


PRECISION  SSEED, SDSEED, UNSEED,  TLSEED,  TSEED 
SSEED, SDSEED, UNSEED, TLSEED, TSEED 

X  MA  X  M (  5  000 ,23)  ,XMINM(5  00  0,  18)  ,0(5  00,5)  ,  INT(5  00,5)  , 

T  RAN  L  (  200 )  ,  INTR(2  00)  ,SUSLN(  1O0)  ,  INTSN(  100)  ,XMC0EF(23,5)  , 

P  S  (  5  )  ,  P  T (  5  )  , S  L  E  N  D  A (  5  )  ,  T  L  E  N  0  A (  5  )  ,  ISLIFE,  NURUN,  NUSPAN,  NUROW, 
ISMEAN, ISDMI N , NDMIN , NDMAX , NT IME , NFMIN , ITMEAN, ITRMIN 
ISPAN, 1RUN , NUQ 1 .NFRAME, NUR0W2 


C 

DIMENSION  S  P  AN  L ( 5  ) 


c 

C  READ  ALL 
C 

READ ( 5 
RE AD ( 5 
READ ( 5 
RE AD ( 5 
READ ( 5 
READ ( 5 
READ ( 5 
READ ( 5 


INPUT  DATA 

101  ) NURUN , S  L I FE , NFRAME ,  I  CHECK 

102  ) NUSPAN  .WIDTH,  RATI  0 

103)  ( S  P  A  N  L  (  I  )  ,  1  =  1  .NUSPAN) 

111) NR0WM1 , NR0WM2 

111) NROWS 1  ,  N ROWS  2 
111) N  R  OWC 1  , N  R  OWC  2 
111) NROWA 1 , NR0WA2 

104 )  NUROW, NUR0W2 


DO  11  1=1, NUROW 

R  E  A  0 (  5  ,  105)  ( X  MC  0  E  F (  I  ,J)  ,J=1  ,  NUSPAN  ) 
11  CONTINUE 


READ ( 5 , 106 ) SUMEAN .SUSTD, AREA, SSEED 
RE  AD  (  5  ,  107)S ME AN.  SDMIN, SDSEED 
READ(5,108)DMIN, DMA X.UNTI ME, FMIN, UNSEED 


1  09 

1  1  0 

1  1  1 

1  1  2 

1  1  3 

i  i  a 

1  1  5 

1  1  6 

1  1  7 

1  1  6 

1  1  9 

1  20 

1  2  1 

122 

123 

1  24 

1  25 

1  26 

127 

1  2  6 

129 

1  30 

1  3  1 

132 

133 

134 

1  35 

136 

137 

1  36 

1  39 

1  40 

1  4  1 

142 

1  4  3 

1  44 

145 

1  46 

1  47 

1  4  8 

1  49 

1  50 

1  5  1 

1  52 

153 

1  54 

1  55 

156 

1  57 

1  56 

1  59 

1  60 

1  6  1 

1  62 

1  63 

1  64 

1  65 

1  6  6 

1  67 

1  66 

1  69 

1  70 

1  7  1 

1  72 

1  73 

1  74 

1  75 

1  76 

177 

1  76 

1  79 

1  60 

1  6  1 

1  62 

1  63 

1  64 

1  85 

1  6  6 

1  67 

1  6  S 

1  89 

1  90 

1  9  1 

192 

1  9  3 

1  94 

1  9  5 

1  9  6 

1  9  7 

1  9  8 

199 

200 

20  1 

202 

203 

204 

205 

206 

207 

208 

209 

2  1  O 

2  1  1 

2  1  2 

2  1  3 

2  1  4 

2  1  5 

2  1  6 
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R  E  A  D  (  5  ,  1  09  )  MR  ,  MRS  TD  ,  XMO , OS  T  D  ,  TLSEED 
R  E  A  D  (  5  ,  1  1 O ) TMEAN ,  TRMIN.TRDUR,  TSEED 
C 

C  WRITE  ALL  INPUT  DATA 

C 

WRITE ( 6 . 1001 ) 

WR  I TE (  6  ,  1 002  ) NFRAME 
WR  I TE {  6 .  1 003 ) 

WR I T  E  (  6  ,  1004  )NURUN , S  L  IFE 
WR  I  TE (  6 ,  1005  ) 

WRITEI6.  1006  )NUSPAN, WIDTH,  RATIO 
WRITE ( 6  .  1007 )  (SPANL (  I  )  .  I  *  1  , NUSPAN ) 

WR  I TE (  6  ,  1008 ) 

WRITEI6. 1020) NROWM 1 , NR0WM2 
WR  I T  E (  6 ,  1021  ) NROWS 1  , N ROWS  2 
WR  I  TE (  6  .  106  1  ) NR OWC 1  .  NR0WC2 
WRITE16.  1022) NROWA 1  ,  NR0WA2 
WR I TE ( 6 , 1 009 ) NUROW, NUR0W2 
C 

DO  12  1*1, NUROW 

WRITE(  6,  1  O 1 O )  ( XMCOEF (  I  ,J)  ,J*  1  .NUSPAN) 

12  CONTINUE 
C 

WR  I  TE (  6 ,  1 O 1  1  ) 

WR  1  TE  (  6  ,  10  12 ) SUMEAN , SUSTO , AREA , SSEED 
WRITE! 6. 1013 )SMEAN , SDMIN , SDSEED 
WRITEI6.  1014) 

WRITE ( 6 , 1015 )DMIN, OMAX , UNTIME , FMIN , UNSEEO 
WR  I T  E (  6  ,  1016) 

WR  I  TE (  6  .  1017JMR, MRS  T  0 ,  XMO, OS TD. TLSEED 
WRITE16,  1016JT ME AN,  TRMIN.TRDUR,  TSEED 
WR  I  T  E {  6 ,  1019) 

C 

1 F ( NURUN  LT  .  2  )  CALL  EX  I  T 
I F ( NUSP AN  .  LT  .  2 ) CAL L  EXIT 
1 F (  I  CHECK  .  EO  .  1  ) CAL L  EXIT 
C 

C  CONVERT  THE  DESIGN  LIFETIME  INTO  THE  PREVIOUSLY  DEFINED  UNITS 
C 

ISLIFE*SLIFE*365  *24 . /TROUR^O . 5 

C 

C  DETERMINE  THE  VARIANCE  OF  SUSTAINED  EUDL  FOR  EACH  SPAN 
C 

WR  I  TE (  S ,  103  1  ) 

SUVAR*SUSTD**2 
SICS2*SUVAR/(0.03+1  .O/AREA) 

S 1 G2  * O  03*S I GS2 
WIDTH  *WI DTH*  2  . 

SIGS2*SIGS2*2  .2/WIDTH 
SOME AN  *SUMEAN*SUMEAN 
C 

DO  2 1  ISPAN* 1 , NUSPAN 

VARS  U  *  S I C  2  +  S I GS2/SPANL  (  ISPAN) 

C 

C  DETERMINE  THE  PARAMETERS  FOR  THE  GAMMA  DISTRIBUTION  OF  THE 
C  SUSTAINED  EUOL 
C 

PS  (  ISPAN)  S SOME AN/ VARSU 
S  L  END A (  ISPAN)  *SUMEAN/VARSU 
C 

C  WRITE  THE  STANDARD  DEVIATION  OF  THE  SUSTAINED  EUDL 
C 

STDSU  *SORT ( VARSU ) 

WRITEI6, 1032)ISPAN, STDSU 
C 

21  CONTINUE 
C 

C  CONVERT  THE  MEAN  AND  MINIMUM  SUSTAINED  LOAD  DURATION  INTO  THE 
C  PREVIOUSLY  DEFINED  UNITS 
C 

I SMEAN : S  ME  A  N  *  3  65  .*24  /TRDUR  +  O.  5 
1SDMIN=SDMIN*365 .*24 . /TRDUR+O . 5 
C 

C  CONVERT  THE  DATA  ON  UNOCCUPIED  PERIODS 
C 

NDMI N : DM] N*24 ./TRDUR+O. 5 
N  D  MA  X  *  D  MA  X  *  2  4  .  /TRDUR  +  O  .  5 
NT1ME  =  UNTIME*365  . *24  .  /TRDUR  +  O  .  5 
NFMI N  s  FMI N*365  .  *24  ,  /TRDUR  +  O  .  5 
C 

C  DETERMINE  THE  MEAN  AND  VARIANCE  OF  TRANSIENT  EUDL  FOR  EACH  SPAN 
C 

WR I TE ( 6 , 1 033 ) 

XMOR*  XMO*FLOAT ( MR ) 

XRSTD2* FLOAT ( MRS  T  D ) **2 

ORSTD : 2 . 2* ( OSTD*OSTD*FLOAT ( MR ) ♦  X  MO  *  X  MO  *XRSTD2*  XMO  R  *  *  2 ) 

DO  22  ISP AN  *  1  .  NUSPAN 
AREA*SPANL( ISPAN)* WIDTH 
I F ( AREA . LT . 30 . 0 ) GO  TO  31 
X  *  ( AREA  -  14.41/0.56 
B  L  END A  *  SORT ( X ) 

GO  TO  32 

31  BLENDAsO. 173*AREA 

32  TRMEAN  *  BL  ENDA*XMQR/ AREA 
TRVAR=BLENDA*ORSTD/( AREA* AREA) 

C 

C  DETERMINE  THE  PARAMETERS  FOR  THE  GAMMA  DISTRIBUTION  OF  TRANSIENT 
C  EUDL 
C 

P  T (  ISPAN)  =TRMEAN*TRMEAN/TRVAR 
TLENDAl  ISPAN)  * T RME AN / T R V A R 
C 

C  WRITE  THE  MEAN  AND  STANDARD  DEVIATION  OF  THE  TRANSIENT  LOAD 
C 

TRSTD  *SORT ( TRVAR  ) 

WR  I T  E (  6 ,  1034  )  ISPAN,  TRMEAN , TRSTD 
C 


2  1  7 

2  1  8 

2  1  9 

220 

22  1 

222 

223 

224 

225 

226 

227 

228 

229 

230 

23  1 

232 

233 

234 

235 

236 

237 

238 

239 

240 

24  1 

242 

243 

244 

245 

24  6 

247 

243 

24  9 

250 

25  1 

252 

253 

254 

255 

256 

257 

258 

25  9 

260 

26  1 

2  62 

2  63 

264 

2  65 

2  6  6 

2  67 

268 

269 

270 

27  1 

272 

273 

274 

275 

27  6 

277 

278 

279 

280 

28  1 

282 

283 

284 

285 

28  6 

287 

288 

269 

290 

29  1 

29  2 

293 

29  4 

295 

29  6 

297 

296 

29  9 

300 

30  1 

302 

303 

304 

305 

306 

307 

308 

309 

3  1  O 

3  1  1 

3  1  2 

3  1  3 

3  1  4 

3  1  5 

3  1  6 

3  1  7 

3  1  8 

3  1  9 

320 

32  1 

322 

323 

324 
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22  CONTINUE 
C 

C  CONVERT  THE  TIME  BETWEEN  OCCURRENCES  OF  TRANSIENT  LOADS  INTO  THE 
C  PREVIOUSLY  DEFINED  UNITS 
C 

I TMEAN  *  TMEAN*24  .  /TRDUR  +  O .  5 
I TRMI N  t  TRMI N  *  2  4  .  /TRDUR  +  O .  5 
C 

WRITE (  6 .  105 1  ) 

WR I  T  E  (  6  ,  1052) NR0WM1  , NR0WM2 
WRITE(6,  105  3 ) NROWS 1  ,  N ROWS  2 
WR I TE (  6  .  1062) NROWC 1  . NR0WC2 
WR I TE ( 6 , 1054 ) NROWA 1 , NR0WA2 
C 


RETURN 


FORMAT  STATEMENTS 


1  O  1 
102 

103 

104 
1  05 
106 
1  07 
1  08 
1  09 
1  1  O 
1  1  1 

1  OO  1 


1  003 
1  004 


1  005 
1  006 


FOR  MAT(I10. MO. 2,2110) 

FORMAT  (  I  1 0 , 2F 10 . 2 ) 

FORMAT ( 5F 1 O . 2 ) 

FORMAT (2110) 

FORMAT ( 5F 1 O . 2 ) 

FORMAT (3F 10. 2,010. 2) 

FORMAT (2F 10.2, DIO. 2) 

FORMAT (4F 10. 2, DIO. 2) 

FORMAT ( 2  I  10,2F10.2,D10.2) 

FORMAT (3F 10. 2, 010. 2) 

FORMAT (2110) 

FORMAT ( '1',////10X, '****••*•***********•*************»**•«*****« 
1***'  / 1 0  X  ,  '  MONTE-CARLO  STUDY  OF  PATTERN  L 0 A 0  I  N C S  (  L I  V E  LOAD)' 

1  / 1 o  X ,  '.«**.«**•«*«*••«****««*»*«******•***********«****'  ) 

FORMAT ( ////5X ECHOCHECK  OF  INPUT  DATA  -  FRAME', 14 
1  /5X ,'**•••**•**«******•*****'  ) 

FORMAT!  //5X,  'GENERAL  DATA'/5X,  ' . ') 

FORMAT (  1  OX  ,' NUMEER  OF  RUNS 
1  /10X, 'DESIGN  LIFE  OF  THE  STRUCTURE 

FOR  MAT ( / / 5  X  ,  'FRAME  PR0PERTIES'/5X,  ' . 

FORMAT (  1  OX  ,' NUMBER  OF  SPANS 
1  /  1  OX ,  '  TR 1 BUTAR Y  WIDTH 

1  /lOX, 'RATIO  OF  COLUMN  TO  BEAM  STIFFNESS* ‘ 

SPAN  LENGTH  IN  ME T E R S  '  /  1  5 X  ,  5 F 1 O .  2 ) 

LOAD  EFFECT  INFLUENCE  COEFFICIENT  MATRIX: 


.110 

.  F  1  O  .  1  , 


TOTAL  NUMBER  OF  ROWS  =',I4 

ROWS  DO  NOT  REQUIRE  LIFETIME  MINIMUM  VALUES) 


1007  F  0  RMA  T (  / 1 0  X 

1008  FORMATI/IOX 

1  /  1  OX 

1009  F  ORMAT (  / 1  OX 
1 / 1  OX ,  '  ( THE  LAST '  ,  I  3 
1  '/) 

1010  FORMAT ( 1 5X , 5F 1 O . 4 ) 

1011  F0RMAT(////5X,  'SUSTAINED  L0AD'/5X,' . 

1012  FORMAT (  1  OX  MEAN  SURVEY  LOAD 

1 R  '  / 1  OX  ,  ' S T ANDARD  DEVIATION  OF  SURVEY  LOAD 

/ 1 O  X ,  'THE  ASSOCIATED  AREA 

'INITIAL  MAGNITUDE  SEED  NUMBER 
'MEAN  LOAD  DURATION 
/  1  OX ,  ' Ml N I  MUM  LOAD  DURATION 
/  1  OX  ,  '  I  N I T I AL  OURATION  SEED  NUMBER 

’UNOCCUPIED  PERI0D'/5X, ' . 

'MINIMUM  OURATION 
'MAXIMUM  DURATION 
'TIME  BETWEEN  TWO  PERIODS 
'MINIMUM  TIME  BETWEEN  TWO  PERIODS 
'INITIAL  SEED  NUMBER 

'TRANSIENT  L  0 A  D  '  / 5  X ,  ' . 

'MEAN  NUMBER  OF  LOADS  PER  CELL  * 

'STANDARD  DEVIATION  OF  LOADS  PER  CELL* 
'MEAN  WEIGHT  OF  SINGLE  LOAD  * 

'STANDARD  DEVIATION  OF  SINGLE  LOAD  * 
'INITIAL  MAGNITUDE  SEED  NUMBER  * 

'MEAN  TIME  BETWEEN  TWO  LOADS  * 

'MINIMUM  TIME  BETWEEN  TWO  LOAOS  = 

'DURATION  OF  ONE  LOAD  * 

'INITIAL  TIME  SEED  NUMBER  * 


. . 

.110 

,  F  1  O  .  2  .  4 
,  F  1  O  .  2  ) 


YEARS ' ) 


METERS ‘ 


) 


1  R  ' 

1  / 1  OX  , 

10  13  FORMAT (  1  OX , 

1 
1 

1014  FORMAT ( //5X , 

1015  FORMAT (  1  OX  , 

1  / 1  OX  , 

1  / 1  OX  , 

1  / 1  OX . 

1  / 1  OX  ,  ‘ 

1016  FORMAT ( //5X , 
10  17  FORMAT (  1  OX , 

1  /  1  OX  , 

1  / 1  OX  . 

1  / 10X . 

1  /  1  OX  . 

10  18  FORMAT (  1  OX  , 

1  / 1  OX  , 

1  / 10X , 

1  / 10X . 


) 


'  ,  F  1  O  .  3 ,  '  KN/SO  METE 
' . F 10 . 3 , '  KN/SO  METE 
' , F 10 . 2 , '  SO. ME TER' 
',018.10) 

'  ,  F  1  O  .  2  ,  ' 

'  ,  F  1  O  .  2  ,  ' 

',018.10) 


YEARS ' 
YEARS ' 


) 


'  ,  F  1  O  .  2  ,  ' 

'  ,  F  1  O  .  2  ,  ' 

'  ,  F  1  O  .  2  ,  ' 

'  ,  F  1  O  .  2  ,  ' 

',018.10) 

-  '  ) 

*  '  ,  I  8 
.  I  8 

,  f  i  o  . : 

,  f  i  o  . : 


DAYS  ' 
DAYS  ' 
YEARS' 
YEARS  4 


,  D 1 6  .  10) 

,  F  1  O  .  2  ,  ' 

,  F  1  O  .  2  ,  ' 

,  F  1  O  .  2  ,  ' 

,018. 10) 


KN  ' 

KN  ' 

DAYS  ' 
DAYS  ' 
HOURS  ' 


10  19 

FORMAT (  / / 5  X  . 

'  *****END 

OF 

INPUT 

DATA*  *  *  *  *  '  ) 

1020 

FORMAT (  / 1  OX  , 

'  ROW 

'.13, 

'  TO 

',13, 

'  CORRESPONDS 

TO 

BEAM 

MOMENTS ' ) 

102  1 

FORMAT (  1  OX  ,  ' 

ROW' 

.13,' 

TO  ' 

,13,' 

CORRESPONDS 

TO 

SHEARS 

'  ) 

1  OG  1 

FORMAT (  1  OX  ,  ' 

ROW  ' 

.13,' 

TO  ' 

.13,' 

CORRESPONDS 

TO 

COLUMN 

MOMENTS ' 

1022 

FORMAT (  1  OX  .  ' 

ROW' 

.13.' 

TO  ' 

,13,' 

CORRESPONDS 

TO 

AXIAL 

FORCES  '  ) 

1C31 

FORMAT (  ///  t  OX .  ' 

SPAN 

STANDARD  DEVIATION 

OF 

SUSTAINED  EUDL 

1  /  3  0  X  ,  ' 

- -KN/SO , METER 

...  j 

1  032 

FORMAT ( 5X . 18 

,  1  7X 

,  F  1  O  . 

3  ) 

1033  FORMAT (// 1 0X ,  '  SPAN  MEAN  OF  TRANSIENT  EUDL 

1  /26X ,' KN/SO . METER  K N / S 0 . ME T E R ' ) 

1034  FORMAT ( 5K ,  I  8  ,  F20 . 3  ,  F23  .  3 ) 

1051  FORMAT! ' 1 ' ,/////////////// 50X, '*************»****•**** 

1  / / 5  O  X  ,  '  RESULTS  OF  ANALYSIS' 

1  //50X, '»**«*«*****************') 

1052  FORMAT ( ////37X ,' FROM  LOCATION 

1METER/METER -WIDTH ' ) 

1053  FORMAT ( ////37X ,' FROM  LOCATION 

1W1 DTH ' ) 

1062  FORMAT ( ////37X ,' FROM  LOCATION 

IN- ME  TER/ METER-WIDTH'  ) 

1054  FORMAT ( ////37X ,' FROM  LOCATION 

1  ME  TER-W10TH '  ) 

END 


13,' 

TO  ' 

.13, 

13.' 

TO  ' 

.13, 

13.' 

TO  ' 

.13, 

13.' 

TO  ' 

.13, 

*  *  *  * 

*  *  *  * 

» 

STO  DEVIATION' 


BEAM  MOMENTS  IN  KN- 


SHEARS  IN  KN/METER- 


COLUMN  MOMENTS  IN  K 


AXIAL  FORCES  IN  KN/ 


SUBROUTINE  SUSTAN 


THIS  SUBROUTINE  SELECTS  THE  NUMBER  AND  MAGNITUDES  OF  SUSTAINED  LOADS 
OCCURRING  IN  ONE  SPAN  IN  THE  LIFETIME  OF  THE  STRUCTURE.  THE  DURATION 


. 


325 

326 

327 

32  8 

329 

330 

33  1 

332 

333 

334 

335 

336 

337 

338 

339 

340 

34  1 

342 

343 

344 

34  5 

34  6 

347 

34  8 

349 

350 

35  1 

352 

353 

354 

355 

35  6 

357 

358 

359 

3  60 

36  1 

362 

3  63 

364 

3  65 

36  6 

3  67 

368 

3  69 

370 

37  1 

372 

373 

374 

375 

376 

377 

37  8 

379 

380 

38  1 

362 

363 

364 

385 

38  6 

367 

38  8 

389 

390 

39  1 

392 

39  3 

39  4 

39  5 

39  6 

39  7 

39  8 

39  9 

400 

40  1 

402 

403 

404 

405 

406 

407 

406 

409 

4  1  0 

4  1  1 

4  1  2 

4  1  3 

4  1  4 

4  1  5 

4  1  6 

4  1  7 

4  1  8 

4  1  9 

420 

42  1 

422 

423 

424 

425 

426 

427 

428 

429 

4  30 

43  1 

432 
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c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


1  3 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


22 

c 


c 

c 

c 


AND  NUMBER  OF  UNOCCUPIED  PERIODS  ARE  ALSO  SELECTED.  LOAD  MAGNITUDES 
( KN/METER/METER -WI DTH ) AND  THE  TIME  AT  THE  BEGINNING  OF  EACH  DIFFERENT 
LOAD  ARE  STORED  IN  'SUSLN'  AND  ' INTSN'  RESPECTIVELY. 

SUBROUTINE  SUSTAN 


DOUBLE  PRECISION  S S E E D , S D S E E D , U NS E E 0 , T L S E E 0 , TS E E D 


COMMON  SSEEO , SDSEED  ,  UNSEED .  TLSEED , TSEED 

COMMON  XMA  X  M ( 5000 , 23 ) . XMI NM( 5000 , 18) .0(500,5) . INT(500,5) , 

1  TRANL  (  200 )  ,  I N  T  R  (  2  O  0 )  ,  S  U  S  L  N  (  100)  ,  I NTSN (  100)  ,XMC0EF(23,5)  , 

1  PS(5)  ,PT(5)  .SLENDAI5)  ,TLENDA(5)  ,  ISL I FE ,  NURUN , NUSPAN ,  NUROW, 

1  ISMEAN, ISDMIN.NDMIN.NDMAX.NT1 ME. NFMIN.IT ME AN.ITRMIN 

COMMON  ISPAN,  IRUN.NU01  , N  FRAME , NUR0W2 


DIMENSION  WK(  1S3)  ,  ISDUR(SO)  4 SUSL  ( 50)  ,R(  1  ) 

ISTOTtO 
NUSUS  e  1 

I SF  x  I SL I FE -  I  SDMI N 
SMEANsFLOAT ( ISMEAN) 

SDMI N  *  FLOAT (  1SDMIN) 


SELECT  SUSTAINED  LOAD  DURATION  FROM  THE  EXPONENTIAL  DISTRIBUTION 
NUR  =  1 

CALL  GGEXN ( SOS EED . SMEAN , NUR .  R  ) 

I  F ( R (  1  )  .  LT  .  SOMI N ) GO  TO  11 


DETERMINE  NUMBER  OF  SUSTAINED  LOADS  IN  THE  DESIGN  LIFETIME 

I  SOUR ( NUSUS  )  :R(  1  ) +0 . 5 
ISTOT  s  I  S  T  0  T  *  ISDUR  ( NUSUS  ) 

IF(  ISTOT.GT  .  I S  F  )  G  0  TO  12 
NUSUS  x  NUSUS-*-  1 
GO  TO  11 

ISDUR(NUSUS)xlSLlFE*  (  ISTOT-ISDUR(NUSL'S)  ) 

SELECT  THE  MAGNITUDE  OF  EACH  SUSTAINED  LOAD  FROM  THE  GAMMA 
D I STR I  BUT  I  ON 


CALL  GGAMR ( SSEED ,  PS (  ISPAN)  . NUSUS, WK, SUSL  ) 
SLEND:SLENDA(  ISPAN) 

00  13  1:1, NUSUS 

SUSL  (  I  )  xSUSL  (  I  )  /S  LEND 

CONTINUE 


NP  x  1 
ISN  xO 
NPSD  x  1 

TIMExFLOAT(NTIME) 

DMINsFLOAT ( NDM1N ) 

DIFFxFLOAT(NDMAX)-DMIN 

00  2  1  ISUS  x  1  , NUSUS 

DETERMINE  NUMBER  OF  UNOCCUPIED  PERIODS  IN  EACH  SUSTAINED  LOAD 
ASSUMING  THAT  AN  UNOCCUPIED  PERIOD  OCCURS  AT  THE  BEGINNING  OF  EACH 
SUSTAINED  LOAD.  THE  MINIMUM  TIME  (NFMIN)  BETWEEN  THE  OCCURRENCE  OF 
THE  LAST  UNOCCUPIED  PERIOD  AND  THE  OCCURRENCE  OF  NEXT  SUSTAINED 
LOAD  MUST  BE  GREATER  THAN  THE  MAXIMUM  DURATION  OF  AN  UNOCCUPIED 
PERIOD  (NDMAX).  NFMIN  MUST  BE  LESS  THAN  THE  MINIMUM  SUSTAINED 
LOAD  DURAT I  ON ( SDMI N  )  . 

ISO:  ISDUR  (  ISUS  ) 

XNUNCxFLOAT(ISD)/TIME-M  .0 
NUNG  =  XNUNG 

ITIMExTIME*  (XNUNG-FLOAT(NUNG)  )  *  O  .  5 
IF( ITIME . LT . NFMINjNUNGsNUNG* 1 

DO  22  J  x  1  ,  NUNG 

ASSIGN  ZERO  LOAD  TO  EACH  UNOCCUPIED  PERIOD  AND  RECORD  THE  TIME  AT 
THE  BEGINNING  OF  EACH  UNOCCUPIED  PERIOD 

ISN: ISN* 1 
SUSLN (  ISN )  sO  .  O 
INTSN (  1 SN  )  x  NP 

SELECT  DURATION  OF  EACH  UNOCCUPIED  PERIOD  FROM  UNIFORM  DISTRIBUTION 
(THE  MINIMUM  DURATION  NDMIN  MUST  BE  GREATER  THAN  OR  EQUAL  TO  THE 
DURATION  OF  A  TRANSIENT  LOAD) 

I NEGxDMIN*DIFF*CGUBFS ( UNSEED )*0 . 5 

RECORD  THE  NON-ZERO  SUSTAINED  LOAD  AND  THE  TIME  AT  ITS  BEGINNING 
IT  IS  ASSUMED  THAT  THE  MINIMUM  DURATION  OF  A  SUSTAINED  LOAD  ISDMIN 
MUST  BE  GREATER  THAN  THE  MAXIMUM  DURATION  OF  AN  UNOCCUPIED  PERIOD 
( NDMAX ) . 

I  SN  x  I SN* 1 

SUSLN ( ISN ) xSUSL ( ISUS ) 

I  NTSN  (  ISN)  :NP* I  NEC 

IT  IS  ASSUMED  THAT  THE  TIME  BETWEEN  TWO  UNOCCUPIED  PERIODS  (NTIME) 
MUST  BE  GREATER  THAN  THE  MAXIMUM  DURATION  OF  AN  UNOCCUPIED  PERIOD 
( NDMAX ) 

NP  x  NP*NT 1  ME 
CONTI NUE 

NPSD*NPSD*  1SD 
NP  x  NPSD 
CONTI NUE 

STORE  THE  LAST  SUSTAINED  LOAD  AT  THE  END  OF  THE  LIFETIME 
I S  N  >  ISN* 1 


< 

433 

434 

435 

436 

437 

436 

439 

440 

44  1 

442 

443 

444 

445 

446 

447 

446 

449 

450 

45  1 

452 

453 

454 

455 

45  6 

457 

456 

459 

4  60 

4  6  1 

4  62 

463 

4  64 

4  65 

4  6  6 

4  67 

4  66 

469 

470 

47  1 

472 

473 

474 

475 

476 

477 

476 

479 

460 

46  1 

482 

483 

464 

465 

48  6 

467 

466 

469 

490 

49  1 

4  9  2 

493 

494 

495 

49  6 

497 

49  6 

499 

500 

50  1 

502 

503 

504 

505 

506 

507 

506 

509 

5  1  O 

5  1  1 

5  1  2 

5  1  3 

5  1  4 

5  1  5 

5  1  6 

5  1  7 

5  1  6 

5  1  9 

520 

52  1 

522 

523 

524 

525 

526 

527 

528 

529 

530 

53  1 

532 

533 

534 

535 

536 

537 

538 

539 

540 
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SUSLN (  ISN )  : SUS  L  ( N  US  US  ) 
I  NTSN (  ISN  )  *  1  SL  I  FE 


RETURN 


C 


END 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  TRANST 


THIS  SUBROUTINE  SELECTS  THE  MAGNITUDES  AND  TIME  BETWEEN  OCCURRENCES 
OF  TRANSIENT  LOADS  IN  ONE  SPAN  IN  THE  LIFETIME  OF  THE  STRUCTURE.  THE 
L  0  A  0  MAGN I TUDES ( KN/METER/METER-W1DTH )  AND  THE  TIME  AT  THE  OCCURRENCE 
OF  TRANSIENT  LOADS  ARE  STORED  IN  'TRANL  '  AND  *  I N  T  R  '  RESPECTIVELY. 

IT  IS  ASSUMED  THAT  THE  OCCURRENCE  OF  TRANSIENT  LOADS  IS  INDEPENDENT 
OF  THE  OCCURRENCE  OF  UNOCCUPIED  P E R I  0 D S  .  T H E R E F 0 R E , A  TRANSIENT  LOAD 
IS  PROBABLE  TO  OCCUR  IN  AN  UNOCCUPIED  PERIOD 


C 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 


1 1 


c 

c 

c 


c 

c 

c 

c 


SUBROUTINE  TRANST 


DOUBLE  PRECISION  SSEED , SDSEED , UNSEED , TLSEED , TSEED 


COMMON 

COMMON 

1 

1 

1 

COMMON 


SSEED  ,  SDSEED  ,  UNSEED ,  TLSEED ,  TSEED 

XMAXMi 5000 , 23 )  , XMINM(5000,  18)  .0(500,5)  ,  INT(5O0,5)  , 

TRANL  ( 200 )  ,  INTR(200)  ,SUSLN(  100)  ,  INTSN(  100)  ,XMC0EF(23,5)  . 

P  S  (  5  )  ,  P  T ( 5  )  , S  L  E  N  D  A ( 5  )  ,  T  L  E  N  D  A ( 5  )  ,  I  S L  I  FE .  NURUN , NUSPAN . NUROW, 
IS ME AN, ISDM1N , NDMI N , NDMAX , NTIME , NFMIN, IT ME AN, ITRMIN 
ISPAN,  1  RUN , NUO 1  , N  F  R  A  ME ,  NUR0W2 


DIMENSION  WK ( 600) ,  R(1) 


ASSUME  THE  FIRST  TRANSIENT  LOAD  OCCURS  AT  THE  BEGINNING 


I  NTR (  1  )  *  1 


I  TOT  *  1 
NUTR  *  2 

ISFs ISL IFE-ITRMIN 
TMEAN:FLOAT( ITMEAN) 
TRMINsFLOAT ( ITRMIN ) 


SELECT  THE  TIME  BETWEEN  OCCURRENCES  OF  TRANSIENT  LOADS  FROM 
EXPONENTIAL  DISTRIBUTION.  THE  MINIMUM  T I  ME (  I T R M I N  ) MU S T  BE  GREATER 
THAN  OR  EQUAL  TO  THE  DURATION  OF  ONE  TRANSIENT  L0AD(1) 


NUR*  1 

CALL  GGEXN ( TSEED . TMEAN , NUR , R ) 

I F ( R ( 1 ) . LT . TRMI N ) GO  TO  11 

DETERMINE  NUMBER  OF  TRANSIENT  LOADS  IN  THE  LIFETIME 

I  NT  I  ME  s  R (  1  ) *0 . 5 

ITOT  * ITOT+INTIME 

I F (  I  TOT  .  CT  .  I SF  )  GO  TO  12 

RECORD  THE  TIME  AT  EACH  TRANSIENT  LOAD  OCCURRENCE 
I  NTR  l  NUTR  )  *  I  TOT 
NUTR*NUTR+1 


C 

C  CONTROL  THAT  THE  NUMBER  OF  TRANSIENT  LOADS  WILL  NOT  EXCEED  THE 
C  STORAGE  SPACE 
C 

I F ( NUTR . EQ . 200 ) CO  TO  12 
C 

CO  TO  11 


c 

c 

c 

c 


c 

c 

C 


1  3 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THE  LAST  TRANSIENT  LOAD  WHICH  IS  ASSIGNED  ZERO  MAGNITUDE  OCCURS  AT 
THE  END  OF  THE  LIFETIME 

I NTR ( NUTR )  *  ISL  I  FE 
TRANL ( NUTR ) *0  O 

SELECT  THE  MAGNITUDE  OF  EACH  TRANSIENT  LOAD  FROM  GAMMA  DISTRIBUTION 
I TR: NUTR* 1 

CALL  GG AMR ( TLSEED , PT ( ISPAN) , ITR.WK, TRANL) 

TLEND = TLENDA ( ISPAN ) 

DC  13  I  s  1  ,  I TR 

TRANL  (  I  )  s  TRANL (  I  )/TLEND 

CONT I NUE 

RETURN 

END 


*s«s*s**s»*ss***ss 
SUBROUTINE  COMBND 


THIS  SUBROUTINE  COMBINES  THE  SUSTAINED  LOADS  AND  TRANSIENT  LOADS. 
THE  LOAD  MAGNITUDES  AND  THE  TIME  AT  THE  BEGINNING  OF  EACH  DIFFERENT 


•  ■ 

54  1 

542 

543 

5  44 

545 

54  6 

547 

546 

5  4  S 

550 

5  5  1 

552 

553 

554 

555 

556 

557 

556 

559 

5  60 

5  6  1 

562 

5  63 

5  64 

565 

5  6  6 

5  67 

566 

5  63 

570 

57  1 

572 

573 

574 

575 

57  6 

577 

576 

579 

5  60 

58  1 

582 

583 

584 

585 

58  6 

587 

588 

589 

590 

59  1 

592 

59  3 

59  4 

595 

59  6 

597 

59  8 

59  9 

600 

60  1 

602 

6  C  3 

604 

605 

606 

607 

608 

609 

6  1  0 

6  1  1 

6  1  2 

6  1  3 

6  1  4 

6  1  5 

6  1  6 

6  1  7 

6  1  6 

6  1  9 

620 

62  1 

622 

623 

624 

625 

62  6 

627 

628 

629 

630 

63  1 

632 

633 

634 

635 

63  6 

637 

638 

639 

640 

64  1 

642 

643 

644 

645 

64  6 

647 

648 
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c 

c 

c 

c 


c 


c 

1  o 
c 


c 

c 

1  1 


c 

c 

c 

c 

c 

c 


c 

c 

1  3 


c 

c 

1  4 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 
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COMBINED  LOAD  ARE  STORED  IN  '0'  AND  '  I N  T '  RESPECTIVELY  . 
SUBROUTINE  COMBND 

OOUBLE  PRECISION  SSEED , SDSEEO , UNSEED , TLSEED , TSEED 


COMMON  SSEED , SDSEED , UNSEED , TLSEED , TSEED 

COMMON  XMAXM ( 5000 , 23  )  , XMI NM( 5  000 ,  18). 0(5  00, 5),  I  NT ( 500 , 5  )  , 

1  TRANL(200)  ,  I N  T  R (  2  00  )  , SUS  L  N (  100)  ,  I  N  T  S  N  (  100)  ,XMC0EF(23,5)  , 

1  PS(5)  ,  PT ( 5  )  , SLENDA ( 5  )  ,TLENDA(5)  ,  ISLIFE,NURUN,NUSPAN,NUROW, 

1  ISMEAN, 1SDMIN.NDM1N, NDMAX ,NTIME,NFMIN, ITMEAN, ITRMIN 

COMMON  ISPAN,  I  RUN  ,  NUO 1  ,  NFRAME  ,  NUR0W2 


N0«  1 
ISN:  1 
I  TRe  1 


CONTINUE 

LNTSNt  I  NTSN (  ISN ) 

LNTR  s  I  NTR (  I TR ) 

I  F (  LNTSN-LNTR)  11,12,13 


0 ( N  0 ,  ISPAN)  s  S  U  S  L  N (  ISN) 
I  NT ( NO ,  I  SPAN )  * LNTSN 
ISNs ISN* 1 
NO  *  NO*  1 
CO  TO  10 


IT  IS  ASSUMED  IN  THE  FOLLOWING  THAT  THE  MINIMUM  DURATION  OF  A 
SUSTAINED  LOAD  MUST  BE  GREATER  THAN  OR  EQUAL  TO  TWICE  THE  DURATION 
OF  A  TRANSIENT  LOAD. 

Q ( N  0 ,  ISPAN)  sSUSLN(  ISN)*TRANL(  ITR) 

I  NT ( NO  .  ISPAN )  *  LNTSN 
I  F (  LNTSN  .  EO .  IS L I FE ) CO  TO  14 
I  NTR2  s  LNTR*-  1 

I  TR  s  ITR*-  1 
ISNs ISN* 1 

NO  s  NO*  1 

I F (  I NTR (  I TR  )  .  EO  .  I  NT R2 ) CO  TO  10 
Q( NO ,  ISPAN  )  s  SUS  L  N (  ISN  -  1  ) 

I  NT ( NO  ,  ISPAN )  *  I NTR2 
NO  *  NO*  1 
GO  TO  10 


0 ( N  0 ,  ISPAN )  sSUSLN (  ISN-1  ) *  T  R AN  L (  ITR) 
I  NT ( NO  ,  ISPAN )  *  LNTR 
I NTR2  e  LNTR* 1 
ITRs  ITR*  1 
NO  *  NO* 1 

I F (  I NTR (  I TR )  .  EQ  .  I NTR2  )  GO  TO  10 
I  F (  INTR2  .  EO .  LNTSN  )  ISN*  I S  N  *  1 
Q(  NO  ,  ISPAN)  : S  US  LN (  ISN-1  ) 

I  NT ( NO  ,  ISPAN )  s  I NTR2 
NO  *  NO* 1 
GO  TO  10 


IF (  ISPAN  .  EQ .  1  ) N  U  0 1  ‘NO 

RETURN 

END 


******************* 
SUBROUTINE  EFFECT 
******************* 


THIS  SUBROUTINE  CALCULATES  LOAD  EFFECTS  IN  EACH  SPAN  WHEN  THERE  ARE 
CHANGES  OF  LOADS  IN  ANY  SPAN.  THE  LIFETIME  (ALGEBRAIC)  MAXIMUM  AND 
MINIMUM  VALUES  PER  UNIT  WIDTH  ARE  SELECTED  AND  STORED  IN  ' X  MAX  M ' 

AND  'XMINM'  RESPECTIVELY 


SUBROUTINE  EFFECT 


DOUBLE  PRECISION  S S E E D , S D S E E D , U N S E E D , T L S E E 0 , T S E E D 


COMMON 

COMMON 

1 

1 

1 

COMMON 


SSEED, SDSEED, UNSEED, TLSEED, TSEED 

X  MA  X  M ( 5  000 ,23)  ,XMINM(5000.  18)  ,0(5  00,5)  ,  1NT(5  00,5)  , 

T  R AN  L  (  200  )  ,  INTR(200)  ,SUSLN(  1O0)  ,  INTSN(  100)  ,XMC0EF(23.5)  , 
PS(5)  ,  P  T ( 5  )  ,  S  L  E  N  0  A ( 5  )  ,  T  L  E  N  0  A ( 5 )  ,  ISLIFE.NURUN.NUSPAN.NUROW, 
ISMEAN, ISDMI N , NDMI N , NDMAX , NTIME , NFMIN , ITMEAN, ITRMIN 
ISPAN, IRUN,NUQ1 , NFRAME , NUR0W2 


DIMENSION  X  MA  X ( 2  3 )  ,XMIN(  18)  ,W(5)  ,K0UNT(5) 

DO  12  I  ROW:  1  ,  NUROW 
XMAX  (  I  ROW)  : O  .  O 
CONTINUE 

BECAUSE  I  NT (  1  ,  ISPAN  )  :  1  FOR  I SP AN =  1  ,  NUS P AN ,  THEREFORE  THE  FIRST  XMIN 
AND  XMAX  CAN  BE  ASSUMED  AS  THE  FOLLOWING. 

DO  3  1  I  ROW:  1  , NUROW 
DO  32  ISPAN: 1 , NUSPAN 

X  MA  X (  I R  OW )  *  X  MA  X (  IROW)*XMCOEF(  IROW,  ISPAN)*0(  1  .  ISPAN) 

CONTINUE 

CONTINUE 


NUR0W1 * NUROW- NUR0W2 


, 


649 

E  5  O 

65  1 

652 

653 

654 

655 

656 

657 

656 

659 

6  60 

6  6  1 

662 

6  63 

6  64 

665 

666 

6  67 

6  66 

6  69 

670 

67  1 

672 

673 

674 

675 

676 

677 

676 

679 

660 

68  1 

662 

683 

664 

685 

68  6 

687 

668 

6  89 

690 

69  1 

69  2 

693 

69  4 

695 

69  6 

697 

69  8 

69  9 

700 

70  1 

702 

703 

704 

705 

706 

707 

708 

709 

7  1  O 

7  1  1 

7  1  2 

7  1  3 

7  1  4 

7  1  5 

7  1  6 

7  1  7 

7  1  8 

7  1  9 

720 

72  1 

722 

723 

724 

725 

726 

727 

728 

729 

730 

73  1 

732 

733 

734 

735 

73  6 

737 

738 

739 

740 

74  1 

742 

743 

744 

745 

74  6 

747 

748 

749 

750 

75  1 

752 

753 

754 

755 

756 


83 
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C 


4  1 


1  1 
C 


1  9 


C 

C 

c 


c 

c 

c 

c 


1  6 
1  5 


C 


C 

c 

c 

c 


25 

34 

C 

1  7 
C 


C 

c 

c 

c 


20 

c 


35 

C 

42 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 


I F ( NUR0W1 . EO . 0 ) GO  TO  41 
DO  33  I ROWt 1 , NUR0W1 
XMI N  (  1  ROW)  *  XMAX (  I  ROW) 

CONTINUE 

CONT I NUE 

DO  11  I  SPAN  *  1  . NUSP AN 
KOUNT (  I  SPAN ) *  2 
CONT I NUE 

K*  KOUNT (  1  ) 

I F ( K . CT  NUO 1 ) CO  TO  13 

DETERMINE  THE  EARLIEST  OCCURRENCE  OF  LOAD  CHANCES  AMONG  THE  SPANS 

MIN* INT(K, 1 ) 

DO  14  ISPAN*2 , NUSPAN 
NO  *  KOUNT (  ISPAN ) 

I F (  I  NT ( NO ,  ISPAN)  .  LT  .  M  1  N  )  M  I  N  *  I  NT ( NO ,  ISPAN) 

CONTINUE 

SELECT  THE  APPROPRIATE  LOADINGS  TO  CALCULATE  LOAD  EFFECTS 
USING  THE  LOAD  EFFECT  INFLUENCE  COEFFICIENTS 

DO  15  ISPAN  *  1  . NUSPAN 
NO  =  KOUNT (  ISPAN ) 

I F (  I  NT ( NO .  ISPAN )  .  EO  . Ml N ) GO  TO  16 
W(  ISPAN  ) *0( NO- 1 ,  ISPAN ) 

CO  TO  15 

W(  ISPAN  ):OlNO,  ISPAN ) 

KOUNT (  I  SPAN )  *  KOUNT (  1 S  PAN )  +  1 
CONT I NUE 

DO  17  J  ROW* 1  , NUROW 
XMOMT *0.0 

DO  1 8  ISPAN  *  1  ,  NUSPAN 

XMOMT * XMOMT +XMCOEF (  I  ROW.  I  SPAN) *W(  ISPAN) 

CONT I NUE 

PICK  THE  GREATER  AND  SMALLER  LOAD  EFFECT.  THEN  FINALLY  DETERMINE 
THE  LIFETIME  (ALGEBRAIC)  MAXIMUM  AND  MINIMUM  VALUES. 

I F (  I  ROW .  GT  .  NUR0W1  ) GO  TO  25 
I F (  XMOMT  .  L T  .  XMI N (  I  ROW)  ) CO  TO  34 
I F (  XMOMT  .  CT .  XMAX  (  1  ROW)  ) XMAX (  I  ROW)  *  XMOMT 
CO  TO  17 

XMI N (  1  ROW)  *  XMOMT 
CONTINUE 
CO  TO  19 
CONTINUE 

STORE  THE  LIFETIME  MAXIMUM  AND  MINIMUM  VALUES 

DO  20  1  ROW*  1  , NUROW 

X  MA  X  M (  I  RUN ,  I  ROW)  *  XMAX (  I R OW ) 

CONT I NUE 

I F ( N  U  R  OW 1  E  0  . O  )  C  0  TO  42 
DO  35  I  ROW*  1  , NUR0W1 
XMINMl  IRUN,  IROW)  *XMIN(  I  ROW ) 

CONTI NUE 

RETURN 

END 


SUBROUTINE  OUTPUT 
******************* 


THIS  SUBROUTINE  SORTS  THE  LIFETIME  MAXIMUM  AND  MINIMUM  LOAD 
EFFECTS  AND  COMPUTES  STATISTICAL  PARAMETERS.  FINALLY  WRITES 
THE  RESULTS  AND  THE  SORTED  VALUES 


SUBROUTINE  OUTPUT 


DOUBLE  PRECISION  S S E E D , S D S E E D , U NS E E D , T L S E E 0 , TS E E D 


COMMON 

COMMON 

1 

1 

1 

COMMON 


SSEED.SDSEEO , UNSEED, TLSEED , TSEED 

XMAXM15 000,23) ,XMINM(5000, 16) ,0(500,5) , INT(500,5) , 

T  RAN  L  ( 200 )  .  1NTR(2  00)  ,SUSLN(  100)  ,  INTSN(  lOO)  ,XMC0EF(23,5)  , 

PS  (  5  )  ,  P  T ( 5  )  ,  S  L  E  N  D  A ( 5  )  ,  T  L  END A (  5  )  ,  ISLIFE, NURUN, NUSPAN. NUROW, 
IS ME AN, ISOMIN , NDMIN , NDMAX, NTIME , NFMIN , IT ME AN, ITRMIN 
ISPAN, IRUN, NUO 1 . NFRAME . NUR0W2 


DIMENSION  RMOMT ( 5000 ) , RM( 99 ) 


DOUBLE  PRECISION  STDV 


XNURUN  *  F  L  OAT ( NURUN ) 
XMULT: XNURUN/ 1 OO . 

MU  LT  *  NURUN/ 100 
DIFF*XMULT-FLOAT(MULT) 

INC*  XNURUNtO  .01+0.5 
MA  X*XNURUN*0. 99+0. 5 
NR  *  NURUN+ 1 
NUR0W1 s NUROW- NUR0W2 


DO  1  1  I  ROW*  1  # NUROW 


757 

7  5  ft 

759 

7  60 

76  1 

762 

763 

764 

765 

766 

767 

766 

769 

770 

77  1 

772 

773 

774 

775 

776 

777 

776 

779 

760 

76  1 

762 

763 

784 

765 

76  6 

787 

768 

769 

790 

79  1 

79  2 

79  3 

79  4 

7  9  5 

79  6 

797 

79  8 

79  9 

800 

60  1 

602 

603 

604 

605 

806 

607 

606 

609 

6  1  O 

6  1  1 

8  1  2 

6  1  3 

8  1  4 

6  1  5 

6  1  6 

8  1  7 

8  1  6 

8  1  9 

620 

82  1 

822 

623 

624 

625 

6  2  6 

827 

8  26 

829 

830 

£  3  1 

832 

633 

634 

635 

836 

837 

638 

639 

840 

8  4  1 

842 

843 

644 

845 

646 

847 

648 

64  9 

650 

85  1 

852 

853 

654 

655 

656 

857 

656 

659 

860 

6  6  1 

662 

663 

664 


84 


c 

c 


c 

c 

c 

c 

c 


1  4 

c 


c 

90 


c 


I F (  I  ROW  GT  .  NUR0W1  ) GO  TO  22 

00  1  2  I  RUN  >  1  . NURUN 

RMOMT (  IRUN)  *  X M I NM (  1  RUN ,  I  ROW ) 

CONTINUE 

SORT  THE  MINIMUM  VALUES  BY  ALGEBRAIC  VALUE 
CALL  VSRTA ( RMOMT , NURUN ) 

CALCULATE  THE  ME A N , S T A N D A R D  D E V  I  A T I  0 N ,  C 0 E F F I C I E N T  OF  VARIATION. 
COEFFICIENT  OF  SKEWNESS  AND  COEFFICIENT  OF  KURTOSIS. 

CALL  STAT ( RMOMT , NURUN , RMEAN , STDV , COVAR , COSKEW. COKUR) 

WRITE  THE  RESULTS  AND  THE  SORTED  LIFETIME  MINIMUM  VALUES 

WR  I TE (  6 .  1 02 )  I  ROW, N FRAME 
WR  I TE (  6 .  1 00 1  ) 

WR I TE ( 6 . 103) RMEAN .STDV, COVAR, COSKEW, COKUR 

IF  NUMBER  OF  RUNS  IS  MULTIPLES  OF  1 OO , VALUES  OF  0.01  FRACTILE  TO 
0.99  FRACTILE  WITH  INCREMENTS  OF  0.01  FRACTILE  WILL  BE  PRINTED 

I  F  (  0  1  FF  .  GE  .  O . 009  ) GO  TO  90 

DO  14  1*1,99 

J  s  NR  -  I NC*  I 
RM|I): RMOMT ( J ) 

CONTI NUE 

WR  I  TE (  6 ,  1 003  ) 

WRITE(6,104)(RM (I). 1*1,99) 

WR I TE (  6 ,  1 005  )  • 

WR I T  E (  6 ,  104)  ( RMOMT (  I  )  ,  1:1  .NURUN) 


22  CONTINUE 

DO  13  I  RUN  :  1  ,  NURUN 
RMOMT (  IRUN)  tXMAXM(  IRUN,  IROW) 
13  CONTINUE 


C 

C  SORT  THE  MAXIMUM  VALUES  AND  CALCULATE  THE  STATISTICAL  PARAMETERS 
C  THEN  WRITE  THE  RESULTS  AND  THE  SORTED  VALUES. 

C 

CALL  VSRTA ( RMOMT , NURUN ) 

CALL  S TAT ( RMOMT , NURUN , RMEAN , STDV , COVAR , COSKEW, COKUR ) 

C 

WR I TE ( 6 ,  1 02  )  I  ROW,  NFRAME 
WR  I TE (  6 ,  1002) 

WRITE! 6, 103) RMEAN, STDV, COVAR, COSKEW, COKUR 
C 

I  F  ( D 1 FF  . GE  . O . 009  ) GO  TO  91 
WR  I TE (  6 ,  1 003 ) 

WR  I  T  E  (  6 ,  10  4)  ( RMOMT (  I  )  ,  I  s  INC , MAX ,  I NC ) 

C 

91  WR I TE ( 6 , 1 005 ) 

WRITE! 6, 104 ) ( RMOMT ( I ) , I s 1 , NURUN) 

C 

102  FORMAT ('  1  ',  //5X ,'  AT  LOCATION*  ,  13,  '  (REFERRING  TO  THE  ROW  NUMBER  OF 
1  THE  LOAD  EFFECT  INFLUENCE  COEFFICIENT  MATRIX)  -  FRAME', 14 
1  /5X  ,'•*****•******«*'  ) 


FORMAT I///10* 

'MEAN 

«  '  ,  F  1  0 . 3 

//10X, 

STANDARD  DEVIATION 

:  '  ,  F  1  0 . 3 

//ioy. 

COEFFICIENT 

OF 

VAR  I  AT  I  ON 

:  '  ,  F  1  0 . 3 

// 1  ox , 

COEFFICIENT 

OF 

SKEWNESS 

s  '  .  F  1  0 . 3 

/  /  1  0  X  , 

COEFFICIENT 

OF 

KURTOS I S 

: ' , F 10 . 3// ) 

104  FORMAT  ( 5  (  2X  ,  1  OF  1 3 . 3/  )  ) 

1001  FORMAT ( //5X ,' L IFETIME  (ALGEBRAIC)  MINIMUM  VALUES  PER  UNIT  WIDTH:' 

1  /  5  X  ,  . . . 

1002  FORMAT ( //5X  ,  '  L I FET 1  ME  (ALGEBRAIC)  MAXIMUM  VALUES  PER  UNIT  WIDTH:' 

1  /  5  X  ,  ' . ') 

1003  cORMAT ( 10X, 'VALUES  OF  0.01  FRACTILE  TO  0.99  FRACTILE  WITH  INCREMEN 

ITS  OF  0.01  FRACT I  LE  :  '  / 1  OX  ,  ' . 

1 . 7) 

1005  FORMAT  (/// 1  OX  ,' COMPLETE  DATA'/IOX,' . '/) 


C 


1 1  CONT I NUE 


C 

C  WRITE  THE  FINAL  SEED  NUMBERS 
C 

WR  I  TE (  6 .  105) 

WR I T  E  (  6  ,  106ISSEED 
WR  I  T  E (  6 ,  107JSDSEED 
WR  I TE (  6 ,  106)  UNSEED 
WR  I  T  E (  6 ,  109  )TLSEED 
WR I TE ( 6 . 110) TSEED 
C 


105 

FORMAT ( //SX , 

,  'FINAL 

SEED  NUMBERS 

106 

FORMAT (  / 1  OX , 

,  'SSEED 

*',018.11) 

107 

FORMAT  (  1  OX .  ' 

’  SDSEED 

»'  ,018.  11) 

108 

FORMAT  (  1  OX  ,  ' 

’  UNSEEO 

*  '  ,  D  1  8  .  1  1  ) 

109 

FORMAT ( lOX, 1 

’  TLSEED 

t  '  ,  D  1  6  .  1  1  ) 

1  1  O 

FORMAT (  1  OX  ,  1 

'TSEED 

:  '  ,  D  1  6  .  1  1  ) 

C 

RETURN 

C 

END 


C 

C 

c  **************** 

c 

c  SUBROUTINE  STAT 

C 

C  **************** 

c 

c 


) 


. 


a  65 

8  6  6 

867 

8  68 

•  69 

•  70 

87  1 

872 

873 

•  74 

876 

876 

•  77 

ft  78 

879 

880 

88  1 

862 

683 

•  64 

86  5 

666 

867 

668 

889 

890 

89  1 

69  2 

893 

894 

69  5 

69  6 

897 

69  8 

699 

900 

90  1 

902 

903 

904 

905 

906 

907 

906 

909 

9  1  0 

9  1  1 

9  1  2 

9  1  3 

9  1  4 

9  1  5 

9  1  6 

9  1  7 

8  1  6 

9  1  9 

920 

921 

9  22 

923 

9  24 

825 

9  26 

9  27 

928 

9  29 

9  30 

93  1 

932 

9  33 

9  34 

935 

9  3  6 

9  37 

9  38 

9  39 

940 

9  4  1 

9  42 

9  43 

9  44 

9  45 

9  4  6 

9  47 

948 

949 

9  50 

9  5  1 

952 

953 

954 

955 

95  6 

957 

9  58 

959 

9  60 

9  6  1 

9  62 

9  63 

9  64 

9  65 

9  6  6 

9  67 

9  66 

9  69 

970 

97  1 

972 


85 


c 

c 

c 

c 

c 

c 


c 


1 1 
c 

c 


1  2 
c 


c 


c 

1  5 

c 


1  3 
C 


c 


1  4 
C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  SUBROUTINE  CALCULATES  THE  ME A N , S T A N D A R D  D E V  I  A T I  0 N . C 0 E F F I C I E N T 
OF  VAR  I  AT  I  ON ,  COEFF I C I  ENT  OF  SKEWNESS  AND  COEFFICIENT  OF  KURTOSIS 

SUBROUTINE  STAT(X,N, XMEAN, STDV, COVAR.COSKEW.COKUR) 

0 IMENS ION  X ( N ) 

DOUBLE  PRECISION  VAR, STDV 

CO VAR  *0.0 
COSKEW*  O . O 
COKURtO  .  O 
XN*FLOAT(N) 

SUM1  *  0  . O 
DO  11  I  •  1  ,  N 
SUM  1  *  S  U  M 1 ♦  X  (  I  ) 

CONTINUE 

XMEAN  *SUM1 /XN 

SUM2  *0  . O 
DO  12  I  *  1  . N 
SUM2  SSUM24-X  (  I  )  *  *  2 
CONTI NUE 

VAR: ( SUM2 • XMEAN*SUM1 )/FLOAT(N* 1 ) 

STDV  =  DSORT ( VAR ) 

IF ( XMEAN . EO . 0 . O ) CO  TO  15 
COVARsSTDV/XMEAN 

IFISTDV  EO.O.OlGO  TO  16 

SUM3  *0 . O 
DO  13  I  *  1  ,  N 
S  U  M3  - S  UM3  *  X (  I  ) • *  3 
CONTINUE 

SKEW*  ( S  UM3  ~  3 . 0*XMEAN*SUM2 ) /  X  N  ♦  2  . O*  (  XMEAN  *  *  3  ) 

COSKEW: SKEW/ ( STDV* *3 ) 

SUM4  *0.0 

DO  14  I  *  1  ,  N 

SUM4  : S  U  M4 ♦ X  (  I  )  *  «  4 

CONTINUE 

XKUR* (SUM4-A . 0*XMEAN*SUM3+6.0* XMEAN* XMEAN*SUM2) /  X  N  *  3  .  O  * (XMEAN**4 ) 
COKUR  *  XKUR/ (STDV**4) 

RETURN 


I  MS  L  ROUTINES 


THE  FOLLOWING  IS  A  BRIEF  DESCRIPTION  OF  THE  I  MS L  SUBROUTINES  AND 
FUNCTIONS  WHICH  ARE  USED  IN  THIS  PROGRAM.  THIS  INFORMATION  IS 
AVAILABLE  IN  THE  I  MS L  LIBRARY  REFERENCE  MANUAL  EDITION  7  JAN.,  1979 
THE  1  MS  L  ROUTINES  ARE  STORED  IN  THE  PUBLIC  FILE  *  I  MS L L  I B 


*  *  *  *  * 

I  MS  L  ROUTINE  NAME  -  CGEXN 

*  *  9  *  * 

PURPOSE  -  EXPONENTIAL  RANDOM  DEVIATE  GENERATOR 

USAGE  -  CALL  G C E X N I DS E E D . XM . N , R ) 

ARGUMENTS  DSEED  -  INPUT.  AN  INTEGER  VALUE  IN  THE  EXCLUSIVE  RANGE 
(1,2147483647).  DSEED  IS  REPLACED  BY  A  NEW  DSEED 
TO  BE  USED  IN  SUBSEQUENT  CALLS.  DSEED  MUST  BE 
TYPED  DOUBLE  PRECISION  IN  THE  CALLING  PROGRAM 

XM  -  INPUT  MEAN  VALUE 

N  -  INPUT  NUMBER  OF  DEVIATES  TO  BE  GENERATED 

R  -  OUTPUT  VECTOR  OF  LENGTH  N  CONTAINING  THE 
EXPONENTIAL  DEVIATES 


I  MS  L  ROUTINE  NAME  -  G G  AMR 


PURPOSE  -  ONE  PARAMETER  GAMMA  RANDOM  DEVIATE  G E N E R A T 0 R , AN D 

USABLE  AS  THE  BASIS  FOR  TWO  PARAMETER  GAMMA 
DEVIATE  GENERATOR 

USAGE  -  CALL  G G AMR ( DS E E D , A , N , WK , R ) 

ARGUMENTS  DSEED  >  INPUT.  AN  INTEGER  VALUE  IN  THE  EXCLUSIVE  RANGE 
(1,2147483647).  DSEED  IS  REPLACED  BY  A  NEW  DSEED 
TO  BE  USED  IN  SUBSEQUENT  CALLS.  DSEED  MUST  BE 
TYPED  DOUBLE  PRECISION  IN  THE  CALLING  PROGRAM. 

A  -  INPUT  SHAPE  PARAMETER  FOR  THE  DESIRED  GAMMA 
FUNCTION.  A  MUST  BE  GREATER  THAN  O. 

N  -  INPUT  NUMBER  OF  DEVIATES  TO  BE  GENERATED 


86 


973 

C 

WK 

- 

VECTOR  OF  DIMENSION  3*(N+1)  USED  AS  WORK  AREA 

9  74 

c 

975 

c 

R 

- 

OUTPUT  VECTOR  OF  LENGTH  N  CONTAINING  THE  GAMMA 

9  7  6 

c 

DEV  I ATES 

977 

c 

976 

c 

9  79 

c 

****** 

9  80 

c 

I  MS  L  ROUTINE  NAME 

- 

GGUBFS 

9  6  1 

c 

****** 

9  82 

c 

963 

c 

PURPOSE 

• 

BASIC  UNIFORM  (0,1)  RANDOM  NUMBER  GENERATOR  - 

964 

c 

FUNCTION  FORM  OF  GGUBS 

9  65 

c 

6  8  6 

c 

USAGE 

• 

FUNCTION  GGUBFS ! DSEED ) 

9  67 

c 

9  6  6 

c 

ARGUMENTS  GGUBFS 

- 

RESULTANT  DEVIATE 

969 

c 

9  90 

c 

DSEED 

- 

INPUT.  AN  INTEGER  VALUE  IN  THE  EXCLUSIVE  RANGE 

9  9  1 

c 

(1,2147463647)  DSEED  IS  REPLACED  BY  A  NEW  DSEED 

992 

c 

TO  BE  USED  IN  SUBSEQUENT  CALLS.  DSEED  MUST  BE 

9  93 

c 

TYPED  DOUBLE  PRECISION  IN  THE  CALLING  PROGRAM. 

9  9  4 

c 

995 

c 

9  9  6 

c 

9  9  7 

c 

I  MS  L  ROUTINE  NAME 

- 

VSRTA 

9  9  8 

c 

***** 

9  9  9 

c 

1  OOO 

c 

PURPOSE 

- 

SORTING  OF  ARRAYS  BY  ALGEBRAIC  VALUE 

1  00  1 

c 

1  002 

c 

USAGE 

- 

CALL  VSRTAIA, LA) 

1  003 

c 

1  004 

c 

ARGUMENTS  A 

- 

ON  INPUT,  A  CONTAINS  THE  ARRAY  TO  BE  SORTED 

1  005 

c 

CN  OUTPUT,  A  CONTAINS  THE  SORTED  ARRAY 

1  006 

c 

1007 

c 

LA 

• 

INPUT  VARIABLE  CONTAINING  THE  NUMBER  OF  ELEMENTS 

1006 

c 

IN  THE  ARRAY  TO  BE  SORTED 

1  009 

c 

10  10 

c 

10  11 

c 

10  12 

c 

10  13 

c 

****  END  OF  THIS  PROGRAM  **** 

10  14 

END 

END  OF  FILE 


